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Abstract 

The Advanced GAmma Tracking Array (AGATA) is a European project to develop and operate the next generation y-ray spectrom- 
eter. AGATA is based on the technique of y-ray energy tracking in electrically segmented high-purity germanium crystals. This 
technique requires the accurate determination of the energy, time and position of every interaction as a y ray deposits its energy 
within the detector volume. Reconstruction of the full interaction path results in a detector with very high efficiency and excellent 
spectral response. The realization of y-ray tracking and AGATA is a result of many technical advances. These include the develop- 
ment of encapsulated highly- segmented germanium detectors assembled in a triple cluster detector cryostat, an electronics system 
with fast digital sampling and a data acquisition system to process the data at a high rate. The full characterization of the crystals 
was measured and compared with detector-response simulations. This enabled pulse-shape analysis algorithms, to extract energy, 
time and position, to be employed. In addition, tracking algorithms for event reconstruction were developed. The first phase of 
AGATA is now complete and operational in its first physics campaign. In the future AGATA will be moved between laboratories 
in Europe and operated in a series of campaigns to take advantage of the different beams and facilities available to maximize its 
science output. The paper reviews all the achievements made in the AGATA project including all the necessary infrastructure to 
operate and support the spectrometer. 

Keywords: AGATA, Gamma-ray spectroscopy, Gamma-ray tracking, HPGe detectors, Digital signal processing. 
Pulse-shape and gaimna-ray tracking algorithms. Semiconductor detector performance and simulations 
PACS: 07.50.Qx, 07.85.Nc, 29.30.Kv, 29.40.Gx, 29.40.Wk, 29.85.Ca, 29.85.Fj 
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1. Introduction 

Contemporary nuclear physics research aims at un- 
derstanding the microscopic and mesoscopic features of 
the nuclear many-body system, determined by the ef- 
fective interactions and underlying symmetries. These 
aims are often addressed by studying the nuclear sys- 
tem under extreme values of isospin, mass, angular mo- 
mentum or temperature. In particular many facets of 
the nuclear system can be probed and understood by 
studying nuclei far from stability. With the inception 
of the new generation of Radioactive Ion Beam (RIB) 
facilities, in the case of Europe FAIR (Darmstadt, Ger- 
many), HIE-ISOLDE (CERN, Geneva, Switzerland), 
SPIRAL2 (Caen, France) and SPES (Legnaro, Italy), 
where a much wider range of unstable proton- and 
neutron-rich nuclei will become accessible, a new era 
is being opened for nuclear physics experiments. 

For five decades, high resolution y-ray spectroscopy, 
in particular with the large germanium detector arrays, 
has become a cornerstone in nuclear structure stud- 
ies. Intense research and development efl'orts during the 
1980's and the 1990's both in Europe and in the USA 
led to the construction of efficient An escape-suppressed 
y-ray spectrometers [1, 2]. Although the escape- 
suppression technique significantly improves the peak- 
to-total ratio in y-ray spectra, it limits the solid angle 
occupied by the germanium detectors and therefore the 
efficiency of the y-ray detection system. This detection 
technique culminated with the design and construction 
of ffie EUROBALL [3, 4] and GAMMASPHERE [5] 
spectrometers in Europe and in the USA, respectively. 
This contributed in a significant way to the impressive 
progress made in nuclear structure research since ffien. 

In recent years, new important technical advances, 
namely that of position sensitive Ge crystals and track- 
ing array technology, were developed to cope with 
Doppler effects due to large source velocities and the ex- 
perimental conditions at the future facilities for intense 
radioactive and high-intensity stable ion beams. These 
conditions are expected to be extremely challenging, re- 
quiring unprecedented levels of sensitivity and count- 
rate capabilities. The required performance figures are 
beyond reach with conventional escape-suppressed ar- 
rays. The use of electrically segmented Ge crystals en- 
ables ffie identification of the individual points of inter- 
action of the y rays within the volume of the Ge crys- 
tals as well as the determination of ffie deposited energy 
with high resolution. Besides the highly-segmented Ge 
detectors, the realisation of such an array requires digi- 
tal sampling electronics to extract energy, time, and po- 
sition information from the detectors' output signals us- 



ing pulse-shape analysis methods. The path of the y rays 
in the Ge crystals can then be reconstructed, making use 
of "tracking" algorithms on the position and energy in- 
formation of the individual interactions, and the full en- 
ergy of the original y ray can be determined. With the y- 
ray tracking technique [6, 7], ffie Compton-suppression 
shields become unnecessary resulting in a large gain in 
efficiency while maintaining spectral quality. Further- 
more, the direction of emission of each individual y ray 
can be determined with high precision, which is cru- 
cial for a good Doppler energy correction and hence to 
achieve a good energy resolution even when y rays are 
emitted from a fast moving nucleus, as is the case in 
most nuclear reactions. 

This radically new concept constitutes a dramatic ad- 
vance in y-ray detection that will have wide-ranging 
applications also in medical imaging, astrophysics, nu- 
clear safeguards and radioactive-waste monitoring, as 
well as establish a new level of detection capability for 
nuclear-structure studies. Given the importance of this 
development and its far-reaching implications, a Euro- 
pean collaboration currently consisting of over 40 in- 
stitutions from 12 countries has been established to de- 
velop and construct a European 4n tracking spectrom- 
eter called AGATA (Advanced GAmma Tracking Ar- 
ray). A similar project, GRETINA/GRETA, is also on- 
going in the USA [7, 8]. The development of a track- 
ing spectrometer in Europe is based on progress made 
within many previous projects, e.g. MINIBALL [9, 10], 
MARS [11], some of which were co-ordinated and sup- 
ported by the EU TMR programme (project title: De- 
velopment of Gamma-Ray Tracking Detectors for An 
Gamma-Ray Arrays) [12, 13]. 

AGATA is a mobile instrument that will move be- 
tween major laboratories in Europe take advantage of 
the range of diff'erent beams and equipment at each lab- 
oratory and to optimize the use of beam time at these 
facilities. AGATA will therefore be operated in a series 
of science campaigns at specific European facilities. It 
is now fully operational in its first physics campaign at 
INFN Laboratori Nazionali di Legnaro (LNL) in Italy, 
utilizing the wide range of stable beams available. Sub- 
sequently it will operate at ffie GSI facility in Germany 
and the GANIL laboratory in France and later at new 
radioactive beams facilities such as FAIR, SPIRAL2, 
SPES and HIE-ISOLDE. The spectrometer will be ex- 
panded over time, in phases, towards the full 4n 180 
detector system. 

This paper describes concisely the AGATA spectrom- 
eter and summarizes all the necessary developments that 
have been performed by the AGATA collaboration for 
its design, construction and operation. These develop- 



3 



ments range from advances in Ge detector technology, 
digital data acquisition systems, signal decomposition 
and y-ray interaction reconstruction, and in many areas 
of the infrastructure needed to support and operate such 
a complex device. 

2. Conceptual design 

The conceptual design of AGATA explored the pos- 
sible configurations of a y-ray tracking array and com- 
pared their performance in a consistent way. This pro- 
cess was described extensively in [14]. Here only the 
basic ideas and the most important results will be sum- 
marised. 

It is evident that, in order to maximise the detection 
efficiency of AGATA, the solid angle coverage should 
be maximised. In addition, to minimise the develop- 
ment and maintenance costs, the solid angle should be 
covered with only a few elementary shapes. The passive 
parts of the array should be minimised by using com- 
posite detectors, implying grouping (clustering) more 
crystals within the same cryostat. Moreover, to simplify 
the handling and maintenance of such complex objects, 
the detectors should rely on the encapsulation technique 
originally developed for the EUROBALL Cluster detec- 
tors [15]. An additional requirement in the conceptual 
design of AGATA was to keep a sufficiently large in- 
ner space inside the array in order to host the foreseen 
complementary instrumentation, which often is indis- 
pensable in the physics programme of AGATA. 

The most elegant way to achieve a large solid an- 
gle coverage with a few elementary shapes relies on a 
decomposition of the icosahedron, namely of the pla- 
tonic polyhedron having the largest number of faces. 
Such a decomposition will always result in 12 regular 
pentagons and a variable number of irregular hexagons 
once projected onto the spherical surface. As discussed 
in more detail in [14], the configurations having 120 or 
180 hexagons were soon identified as the most attrac- 
tive ones for AGATA, given the possibility to cover the 
solid angle with a few crystal shapes (2 and 3 for the 
case of 120 and 180 hexagons, respectively) and to eas- 
ily form clusters of crystals (one kind of cluster with 4 
and 3 crystals for the case of 120 and 180 hexagons, 
respectively). The contribution to the overall detection 
efficiency provided by the 12 pentagonal detectors was 
considered too limited to justify the extra costs for de- 
velopment. Therefore, the geometry of the array was 
optimised by minimizing the size of the pentagons. In- 
deed, the pentagonal holes can be utilised for mechan- 
ical support, insertion of complementary detectors and 
for the beam entry and exit pipes. 



The possible configurations for AGATA were eval- 
uated through detailed Monte Carlo simulations of the 
full array. The simulation code for AGATA is based 
on the C++ classes of Geant4 [16], which provide a full 
description of the microscopic interactions of radiation 
with matter, as well as tools to implement the geometry 
of complex detector arrangements and to process and 
extract the relevant information. The Geant4 geometry 
libraries were complemented with a specific class capa- 
ble of describing irregular convex polyhedral shapes as 
the AGATA elementary shapes. 

Establishing the optimal geometrical configuration 
for AGATA was a complex problem where the ingre- 
dients considered went beyond the mere overall perfor- 
mance figures of the array and other factors such as the 
reUability, simplicity, symmetry and cost of the adopted 
solution were taken into account. 

From the results of the Monte Carlo simulations de- 
scribed in [14], it can be concluded that the configura- 
tion with 180 hexagonal crystals (Fig. 1) has better en- 
ergy resolution, full-energy efficiency (Fig. 2) and peak- 
to- total (P/T) ratio than any configuration based on 120 
detectors in a broad range of experimental conditions. 
The difference in performance is particularly evident at 
high y-ray multiplicity, see Fig. 3. The performance fig- 
ures for 1 MeV photons of the final optimised geometry 
with 180 detectors are 82 % solid angle coverage, 43 % 
(28 %) full-energy efliciency and 59 % (43 %) P/T ratio 
at a photon multiplicity My - 1 {My - 30). The fi- 
nal configuration for AGATA was therefore chosen to 
be based on 180 segmented hexagonal crystals. 

With hexagonally shaped crystals the azimuthal seg- 
mentation is quite naturally based on 6 sectors, each of 
them centred on the crystal corners. The pattern of the 
longitudinal segmentation was optimised on the basis 
of detailed electric field simulations. On one hand, the 
eff'ective volume of the segments should be balanced. 
On the other hand, the number of longitudinal segments 
should be large enough to achieve the required posi- 
tion resolution after the pulse-shape analysis process, 
namely 5 mm FWHM as calculated in the Monte Carlo 
simulations. A longitudinal segmentation based on 6 
"rings" was considered the best compromise between 
the required performance and the cost of the associated 
electronics. For further details regarding the crystal seg- 
mentation scheme, see subsection 3.1. 

3. Detectors 

The AGATA detectors are based on encapsulated and 
electrically segmented closed-end coaxial n-type high- 
purity germanium (HPGe) crystals. The crystals have a 



4 



Figure 1 : (Colour online) Computer aided design images of the tiling 
of the sphere (left) and the 180 crystal configuration (right). The 
cryostats and the detector encapsulation are not shown. 




Energy (keV) 

Figure 2: Simulated full-energy efficiency for the 120 crystals and the 
180 crystals (AGATA) arrays as a function of the y-ray energy and at 
multiplicity My = 30, emitted by a point source recoiling along the 7, 
axis, with velocity v/c = 0, and vjc = 0.5. 



tapered hexagonal geometry with an asymmetric shape 
to fit into the An 1 80 detector geometry (see Fig. 1 ). This 
geometry is reahsed with three different shapes, with a 
triplet of crystals arranged in identical triple cryostats, 
such that the full An array has 60 cryostats. With this 
configuration, having a 9 cm thick germanium shell, a 
solid angle coverage of up to 82 % is realised. Its rather 
large inner radius of 22.5 cm (to the endcap face of the 
cryostat) allows the use of most ancillary detectors. 

The AGATA triple cluster (ATC) detector contains 
three 36-fold segmented HPGe crystals. The total en- 
ergy deposited in each crystal is collected in the central 
contact (core) leading to 37 signals per crystal. There- 
fore, the ATC detector contains 111 high-resolution 
spectroscopy channels. A more detailed description of 
the ATC detector is given in [17, 18]. All signal chan- 
nels are equipped with a cold preamplifier stage op- 
erated close to the liquid nitrogen temperature of the 
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Figure 3: Simulated full-energy efiiciency and peak-to-total (P/T) ra- 
tio for the 120 crystals and the 180 crystals (AGATA) array for 1 MeV 
photons as a function of the y-ray multiplicity. The y rays have been 
emitted by a stopped point-like source. 



cryostat. The core preamplifier is characterised by low 
noise and a large dynamic range for energy detection, 
pulse-shape analysis and timing properties. A novel re- 
set technique of the core preamplifier allows for an in- 
creased counting rate capability of the detectors of more 
than 50 kHz, preserving an energy resolution close to 
the nominal one. Moreover, the energy range of the 
ATC detector is substantially extended from 20 MeV 
to 180 MeV, see subsection 3.3. These developments 
are documented in [19-21]. Despite the high electronic 
integration density only small cross-talk contributions, 
typically less than 1 in lO""*, are measurable between 
the segments within a crystal. This cross-talk contribu- 
tion is an expected effect caused by capacitive coupling 
of the signals via the bulk Ge material and can be well 
described within an electronic model of the combined 
crystal and preamplifier assembly [22, 23]. 

3.1. The AGATA crystals 

All detectors are produced by the company Can- 
berra, France. The three types of detectors employed in 
AGATA merely differ in their irregular hexagonal shape 
(see Fig. 4). The different geometries are assigned a 
letter and a colour: A - red, B - green and C - blue. 
A serial number is also assigned to each crystal (AOOl, 
A002, etc.). The crystals have a length of 90±1 mm and 
a diameter of 80^° ] mm at the rear At the front they are 
tapered to a hexagonal shape with a 10° tapering angle. 
The crystal's central hole has a diameter of 10 mm and 
extends to 13 mm from the front end. The 6-fold sector- 
wise segmentation goes through the middle of each flat 
hexagonal side. The 6-fold longitudinal segmentation 



forms rings of 8, 13, 15, 18, 18 and 18 mm in thick- 
ness starting at the hexagonal front face of the crystal 
(see Fig. 4). The thicknesses of the rings have been op- 
timised for a uniform distribution of the y-ray interac- 
tions and optimal pulse-shape sensitivity [11]. The seg- 
ment labelling scheme of the AGATA crystals is shown 
in Fig. 5. 




Figure 4: (Colour online) Drawing of the three AGATA crystal ge- 
ometries. The AGATA triple cluster detector combines the three dif- 
ferent crystal shapes. The side view (lower right) shows the position 
of the segmentation lines. All dimensions are given in mm. 



The weight of a bare AGATA Ge crystal is about 2 kg. 
All crystals are made of «-type HPGe material with an 
impurity concentration specified to be between 0.4 and 
1.8 X 10'" cm"-'. The surfaces of these crystals are very 
delicate and therefore each crystal is encapsulated into 
a hermetically sealed aluminium canister with a 0.8 mm 
wall thickness (see Fig. 5). The encapsulation tech- 
nology was developed for the EUROBALL cluster de- 
tectors [24] and extended to segmented detectors within 
the frame of the MINIBALL project [9]. The distance 
between capsule walls and crystal side faces is 0.4 mm 
to 0.7 mm. The 6x7 connector feed-throughs provide 
access to each of the 36 segmented outer contacts. The 
core contact, which is used for applying the high volt- 
age and to obtain the core energy signal, is isolated with 
ceramic material. Efficient y-ray tracking requires ex- 
treme care with the close packing of the crystals. The 
capsules are accurately mounted in the ATC detector 
with a 0.5 mm spacing between the flat surfaces. 




Figure 5: (Colour online) Segment labelling scheme of the AGATA 
HPGe capsules. Along the crystal axis the external contact is subdi- 
vided into six rings labelled 1 to 6. Each ring is subdivided into six 
sectors labelled a to f. 

3.2. The AGATA cryostats 

The cryostats of the ATC detectors were assem- 
bled and successfully commissioned in a common ef- 
fort by the company CTT Montabaur, together with the 
AGATA collaboration. 

The preamplifiers of segment and core contacts are 
divided in two spatially separated parts. The cold input 
stages of the preamplifiers are operated close to the Ge 
crystals. Cooling and mounting in close proximity to 
the detector is required to optimise noise performance. 
In addition a good electronic shielding between the in- 
put stages is required in order to minimise cross-talk ef- 
fects. The AGATA cryostats employ a separated cooling 
scheme for the encapsulated Ge detector and the cold 
part of the preamplifier electronics. While the Ge de- 
tectors are cooled to 90 K, the FETs (Field Effect Tran- 
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sistors) are operated at temperatures near 130K where 
their noise contribution is minimal. The other adjacent 
parts of the preampHfier electronics contribute less to 
the noise performance and are therefore situated outside 
the vacuum, where they are readily accessible. The elec- 
tric connection between the two parts is made by several 
hundreds of individual thin wires with low thermal con- 
ductivity. 

The thermal isolation is established by a vacuum with 
pressure values below 1 x lO^^mbar. This pressure is 
maintained over long periods by the active getter mate- 
rials built into the cryostat. 

Although each individual FET has only an electric 
power consumption of ~20 mW, the total consumption 
of the 111 FETs in one ATC detector adds up to 2.3 W. 
Together with the enhanced thermal connection by the 
wiring inside the cryostat and the radiative heat absorp- 
tion, a considerable cooling capacity is demanded. The 
Dewar of the triple cryostat contains up to 4.5 1 of liquid 
nitrogen. It has a length of 38 cm and an outer diam- 
eter of 25 cm. A full Dewar is sufficient for about 8 h 
of continuous operation. An electronic measurement of 
the liquid nitrogen filling level, which is based on a ca- 
pacitance measurement between a metallic cylindrical 
tube inside the Dewar and the inner wall of the cryo- 
stat [18], is incorporated in the Dewar. The temperature 
is monitored by two platinum resistance thermometers 
of the type PTIOO, one attached to the copper cooling 
finger close to the Dewar and the other one positioned 
close to the crystals. 

The triple cryostats have a length of 92 cm and a 
weight of 48 kg including the crystals and liquid nitro- 
gen. Very tight tolerances are demanded for the man- 
ufacturing of the cryostat endcaps such that the final 
spacing between the endcap side faces of neighbouring 
triple cryostats is 0.5 mm. Fig. 6 shows 5 ATC detectors 
mounted into the support structure at LNL and demon- 
strates the challenges in the design, assembly and on- 
site installation of such cryostats. 

3.3. Preamplifiers 

The preamplifiers for the AGATA detectors require, 
besides the traditional good energy and timing prop- 
erties, also fast and clean transfer functions to regis- 
ter unperturbed signal traces for pulse-shape analysis. 
In addition, a high count-rate capability was demanded 
in order to exploit fully the high geometrical efficiency. 
New preamplifiers have been developed by the AGATA 
collaboration which fulfill these requirements [19-21]. 
The segment and core signals of the AGATA detectors 
are read out through advanced charge-sensitive resis- 
tive feed-back preamplifiers, employing a new fast-reset 




Figure 6: (Colour online) Photograph of the setup with five AGATA 
triple cluster detectors installed at LNL in Italy. 

technique for dead time and dynamic range optimisation 
as well as an original circuit structure for maximising 
the open-loop gain of the charge-sensing stage. 

The preamplifiers have a cold and a warm part. The 
cold part consists of a low-noise silicon FET, model 
BF862, a 1 .0 pF feedback capacitance and a 1 GO feed- 
back resistance. A dedicated shielding was developed 
for the cold preamplifier board for minimising the inter- 
channel cross talk. The warm part, operated at room 
temperature, is separated by 15 cm cabling inside the 
vacuum from the cold part of the cryostat and comprises 
a low-noise transimpedance amplifier, a pole-zero stage, 
a differential output buffer, and fast-reset circuitry. 

The core preamplifier [21] is characterised by low 
noise and a large dynamic range for energy detection, 
pulse-shape analysis and timing properties. Transient 
signals are not deformed due to the large bandwidth. 
The core preamplifier possesses high count-rate capa- 
bilities and an embedded precision pulser. 

The fast reset is provided by the desaturation cir- 
cuitry, which is capable of detecting saturated signals. 
In such a situation a current source is connected, dis- 
charging the capacitance in the pole-zero network and 
achieving a fast restoration of the output level. 

This innovative time-over-threshold (TOT) technique 
[20] for high-resolution spectroscopy extends the y-ray 
energy range up to about 180MeV. In this technique 
the energy is obtained through a precise determination 
of the reset time, required by the desaturation circuitry, 
which is measured as a time difference. The time mea- 
surement is started when the electronic pulse exceeds 
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the saturation threshold and is stopped when the pulse 
amplitude becomes smaller than this threshold. The 
time difference is strictly related to the pulse height of 
the saturating energy signal and allows an energy de- 
termination with a good energy resolution (FWHM) of 
about 0.2 %. This TOT energy resolution value can 
be compared for energies around lOMeV with values 
of about 0.15 % achieved by the standard pulse height 
mode. Beyond 10 MeV the FWHM values are even 
comparable in both modes. 

A custom programmable high-precision pulser is lo- 
cated on the core-preamplifier board. Its applications 
are: testing, calibration, time alignment and dead-time 
corrections, which are relevant for efficiency measure- 
ments of the detector The pulser is used to inject cal- 
ibration pulses to the core electrode itself as well as to 
all segment electrodes through the detector bulk capac- 
itance. The output signal of the pulser is DC coupled to 
the source pin of the core input FET through a resistor 
divider consisting of a 48.5 resistor and a grounded 
1 .8 Q resistor. Thereafter the signal reaches each of the 
36 detector segments via the capacitive coupling of the 
core to the segments. 

A detailed description of the newly developed seg- 
ment preamplifiers is given in [21]. Three segment 
preamplifier channels are integrated on one printed cir- 
cuit board. The power consumption per segment chan- 
nel is limited to 350 mW allowing the 108 closely 
packed spectroscopic channels to be operated close to 
the vacuum feed-throughs in air 

Differential signal outputs of the 111 spectroscopic 
channels are transmitted through 21 MDR (Mini D Rib- 
bon) high-speed digital data-transmission cables. The 
segment MDR connectors merge the 6 segment signals 
of each sector. The core preamplifier has an individual 
MDR connector for the preamplifier output signal and 
the pulser control signals. Table 3 of ref. [17] sum- 
marises the most relevant specifications and their val- 
ues. 

3.4. AGATA detector specifications 

The individual AGATA detectors have to meet the 
detector specifications given by the collaboration. The 
core energy resolution (FWHM) is specified to be better 
or equal to 2.35 keV (1.35 keV) at 1.33 MeV (122keV) 
and the peak shape FWTM/FWHM (FWTM = full 
width at tenth maximum) is smaller than 2.00. The 
segment FWHM at 1.33 MeV (60keV) is specified to 
be better or equal to 2.30 keV (1.30keV) with a mean 
value of the 36 segments values better or equal to 
2.10keV (1.20 keV). The cross talk between channels 
has to be smaller than 1 x 10 Upon delivery, the 



specifications are verified during a customer acceptance 
test, which is performed by the AGATA collaboration 
at three sites, IKP Cologne, University of Liverpool 
and CEA Saclay. For these measurements single test 
cryostats are used, which are equipped with 37 cold 
input stages. Typically standard analogue commercial 
electronics is used for energy-resolution measurements. 
The cross-talk properties are determined with a 37 chan- 
nel coincidence electronics based on high-speed digital 
sampling electronics. They are extracted from a ^"Co 
measurement after adding the coincident signals of any 
pair of segments as the variation of the 1332.5 keV full- 
energy peak position. The peak shift should not exceed 
0.65 keV. The most relevant parts of the specifications 
are summarised in detail in Table 1 of ref. [17]. 

3.5. Performance of the AGATA triple cluster detector 

The ATC detector is equipped with three single core 
preamplifier boards and 3x12 triple segment preampli- 
fier boards. Energy resolutions at 60keV for all seg- 
ments and at 122 keV for the cores, are measured with 
analogue electronics. At higher energies (1332.5 keV), 
the measurements are performed also with digital elec- 
tronics. The results obtained with a triple cryostat are 
compared in Fig. 7 at low energy and in Fig. 8 at 
high energy with the measured performance of the same 
crystals in a single test cryostat. Average values of the 
energy resolution measurements for the first five ATC 
detectors are summarised in Tables 1 and 2. 

The values obtained for the segments in the triple 
configuration are on average even better than in the sin- 
gle test cryostat. Since these resolutions are dominated 
by electronic noise, it demonstrates the successful de- 
sign and integration of the new AGATA triple cluster de- 
tector. Especially the electronic properties comprising 
the cold and warm parts of the new AGATA preamplifier 
assembly is causing very low noise contributions in the 
triple cryostat despite the high integration density of 1 1 1 
analogue channels. In addition, an improved ground- 
ing was applied as a result of various iterations during 
the project. This reduces unwanted high-frequency and 
noise components and brings the energy resolution of 
the ATC detector at low energy to a value well within 
the specification. 

The cross-talk contributions were investigated by 
analysing the coincident traces over a 7 (,is-long time pe- 
riod using the digital acquisition system on all three de- 
tectors. After identification of the true energy deposi- 
tion in exactly one segment the coincident and simulta- 
neous baseline shifts, which occur in all remaining 107 
nonhit segments, are recorded for these one-fold events. 
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Table 1 : Energy resolution values for the first five triple cluster detectors measured with analogue electronics at IKP Cologne. The core FWHM 
values were measured at 122 keV. The segment's average values and their standard deviations were measured at 60keV. Not all the measurements 
were performed for ATC5. 







Core FWHM [keV] 


Segment average FWHM [keV] 


Detector 


Crystal 


Single 


Triple 


Single 


Triple 






cryostat 


cryostat 


cryostat 


cryostat 




AOOl 


1.34 


1.44 


1.08 ± 0.07 


1.01 ± o.os 


ATCl 


B002 


1.29 


1.41 


1.09 + 0.09 


1.04 ±0.07 




C002 


1.28 


1.21 


1.03 ±0.08 


0.97 ± 0.06 




A003 


1.22 


1.42 


1.14 ±0.08 


1.05 ±0.07 


ATC2 


B003 


1.28 


1.36 


1.06 ± 0.07 


1.00 ±0.07 




COOS 


1.24 


1.49 


1.16 ±0.07 


1.14±0.11 




A002 


1.26 


1.44 


1.03 ± 0.08 


0.93 ±0.11 


ATC3 


BOOS 


1.08 


1.43 


l.OS ± 0.08 


l.OS ± 0.08 




C006 


1.09 


1.42 


LIS ±0.10 


1.14 ±0.11 




AOOS 


1.23 


1.28 


1.04 ±0.09 


1.11 ±0.17 


ATC4 


BOOl 


1.29 


1.27 


1.02 ± 0.07 


1.03 ±0.10 




C003 


1.16 


1.33 


1.00 ± 0.09 


1.11 ±0.38 




A004 


1.27 


1.21 


1.17 ±0.08 




ATC5 


B009 


1.36 


1.S4 


1.11 ±0.07 


l.OS ±0.11 




C004 


1.30 




1.11 ±0.08 





Table 2: Energy resolution values (FWHM) at 1332.5 keV for the first five AGATA triple cluster detectors. The single cryostat measurements were 
performed at IKP Cologne with analogue electronics as was the measurements marked ATC Cologne. The measurements at LNL were performed 
with the ATC detectors mounted on the frame using the AGATA digital electronics and data acquisition system. Average values and their standard 
deviations are given for the segments. Not all the measurements were performed for ATC5 at Cologne. 



Detector 


Crystal 


Core FWHM [keV] 


Segment average FWHM [keV] 


Single ATC at ATC at 
cryostat Cologne LNL 


Single ATC at ATC at 
cryostat Cologne LNL 


ATCl 


AOOl 
B002 
C002 


2.33 2.46 2.S0 
2.27 2.46 2.43 
2.2S 2.33 2.42 


2.09 ±0.16 2.19 ±0.10 2.01 ±0.13 
2.13 ±0.11 2.10 ±0.14 1.99 ±0.09 
2.03 ±0.12 2.11 ±0.12 1.94 ±0.11 


ATC2 


A003 

B003 
COOS 


2.28 2.41 2.S6 

2.23 2.S2 2.42 
2.20 2.21 2.39 


2.10 ± 0.13 2.06 ± 0.08 2.06 ± 0.10 

2.08 ±0.11 2.02 ±0.09 1.94 ±0.08 
2.21 ±0.09 2.21 ±0.08 2.08 ±0.11 


ATC3 


A002 
BOOS 

C006 


2.31 2.40 2.S2 
2.29 2.42 2.49 

2.16 2.27 2.S8 


2.07 ±0.11 2.02 ±0.09 1.98 ± 0.09 
2.09 ±0.14 2.13 ±0.11 2.04 ±0.13 

2.12 ±0.09 2.09 ±0.09 2.13 ±0.1S 


ATC4 


AOOS 
BOOl 

C003 


2.23 2.40 2.19 
2.17 2.S0 2.30 

2.34 2.3S 2.40 


2.03 ±0.10 2.08 ±0.13 1.91 ±0.11 
2.06 ±0.11 2.04 ±0.11 1.91 ±0.11 
2.08 ±0.11 2.08 ±0.09 2.04 ±0.21 


ATC5 


A004 
B009 
C004 


2.31 2.36 2.33 
2.33 2.49 2.63 
2.23 2.26 


2.10 ±0.11 2.04 ±0.12 
2.03 ±0.14 1.96 ±0.14 
2.17 ± 0.10 2.04 ± 0.24 
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Figure 7: (Colour online) Energy resolution values of the crystals 
A003, B003 and C005 for the core signals at Ey = 122 keV and for 
the segment signals at Ey = 60 keV. The measurements were done at 
IKP Cologne using standard analogue electronics. The filled circles 
are the results of measurements performed with the crystals mounted 
in a single test cryostat while the crosses show measurements per- 
formed with the crystals mounted in a triple cluster detector (ATC2). 
The specification limits at 122 keV and 60keV are shown as solid 
lines. 



The correlation between the energy deposition in a sin- 
gle detector segment and the energy shift in all other 
segments is determined over an energy range given by 
y-ray emission fi^om ^"^Co and '■'^Cs sources. The ob- 
served cross talk is within the specifications and cross- 
talk contributions are only observed for segment com- 
binations within the same detector crystal. The regular 
pattern is mainly caused by the diff'erent capacitances 
between the core and segment electrodes. A subset of 
105x105 combinations is shown in Fig 9. 

The method applied to quantify this result is based 
on all possible combinations, shown in Fig. 9 for an 
ATC detector with the 11 Ox 1 11 possible cross-talk ma- 
trix elements within the full triple cryostat. The cross- 
talk contributions between segments of different detec- 
tors are on the lO""* to 10"^ level, which is well within 
the acceptable Umits and which can be disregarded for 
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Figure 8: (Colour online) Energy resolution values of the crystals 
A003, B003 and C005 for the core and segment signals at Ey = 
1332.5 keV. The filled circles and crosses show results of measure- 
ments performed at IKP Cologne using standard analogue electronics 
and with the crystals mounted in a single test cryostat and in a triple 
cryostat (ATC2), respectively. The filled squares show results of mea- 
surements performed at LNL with ATC2 mounted on the frame and 
by using the AGATA digital electronics and data acquisition system. 
The energy specification limit at 1332.5 keV is shown as a solid line. 



a standard operation of the ATC detector. 

The new method to determine precise and absolute 
cross-talk matrix elements was applied to all available 
AGATA triple detectors. The cross-talk pattern was 
measured to be at the 0.1 % level, comparable with the 
values shown in Fig 9. The observed structure can be 
entirely attributed to the capacitive coupling between 
core and segments via the bulk Ge material. This be- 
havior is well reproduced by a linear electronic model. 
This cross-talk contribution is considered as the inher- 
ent cross-talk limit given by the construction of the new 
detectors [23]. The calculated values are closely ap- 
proached by the obtained results. The cross-talk contri- 
butions in all AGATA detectors behaved very similarly. 
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implying that the development of the AGATA cryostats 
and the preamplifier electronics has progressed in such a 
way that the fundamental constraints are accounted for, 
detectable and understood. 




Hit Segment 



Figure 9: (Colour online) Observed relative cross-talk contributions 
of nearly all segment combinations between the three crystals of the 
ATC2 detector are plotted (only 105 channels were operational at the 
time of the measurement). The colour scale gives the relative cross 
talk in units of 10"^. A relative energy shift on the ~10"^ level is 
caused by cross talk for segment combinations within one of the three 
detectors. Cross talk between the three different detector capsules is 
observed to be negligible. The other tested ATC detectors show simi- 
lar results. 



Such knowledge provides the project with the informa- 
tion necessary to enable pulse-shape analysis and y-ray 
tracking. 

The AGATA collaboration has two operational ex- 
perimental characterisation centres, based at the Uni- 
versity of Liverpool in the UK and at CSNSM Orsay 
in France. Three new centres are being commissioned 
at GSl Darmstadt in Germany, at IPHC Strasbourg in 
France and at the University of Salamanca in Spain. 
This would significantly increase the number of AGATA 
detectors that could be experimentally characterized. 

4.1. Liverpool scanning system 

A schematic diagram outlining the University of Liv- 
erpool AGATA detector scanning system is shown in 
Fig. 10 [25]. The figure displays the mechanical con- 
figuration of the system, with the detector vertically 
mounted above the colhmated source assembly. A 
920 MBq '^^Cs source is mounted at the end of a 
1 mm diameter coaxial tungsten collimator of 120 mm 
in length. The collimator is mounted in a lead collar and 
source housing assembly, which shields the system op- 
erators from the mounted source. The whole assembly 
is mounted on a precision Parker x-y positioning table. 
The table is moved in precise computer-controlled steps 
through the use of Pacific Scientific stepper motors and 
two Parker Automation axis indexers. The system has a 
position accuracy of 100 (im and can scan over an area 
of the size 30 cm x 30 cm. 



4. Detector characterisation 

The success of the AGATA project relies on the abil- 
ity to reconstruct the trajectories of y rays scattered 
within the germanium detectors. The deposited en- 
ergy and the location of the photon interactions can 
be extracted from the measurement and the analysis of 
the waveforms that arise on the segments during the 
charge collection. The pulse-shape analysis algorithms 
currently developed use databases of calculated pulse 
shapes. These calculated pulses need to be validated 
with real pulse shapes taken at various points within a 
detector. For the development of tracking algorithms, it 
is also crucial to determine experimentally the interac- 
tion position sensitivity in three dimensions in the whole 
volume of a detector. 

The AGATA collaboration has performed a detailed 
analysis of the response function of the crystals. This 
work included developing a theoretical basis data set, 
which describes the detector response function and then 
validating this against the equivalent experimental data. 




— Support stmcture 
Dewar 
LVDS cables 
Detector cryostat 
Scatter collimators 
Scintillators 
Steel plate 
Pb collar 

Injection collimator 

Source housing 
Parker scanning table 

Scintillator 
pre- amp box 
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Figure 10: (Colour online) Schematic diagram of the University of 
Liverpool scanning table assembly. 



A steel frame is constructed around the scanning ta- 
ble. From this, a steel plate is suspended by threaded 
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rods, on which the scatter colUmators and scatter detec- 
tors can be supported. The rods enable the plate, and 
hence the collimators, to be moved in the z-direction. 
The steel plate has a square cut from its centre and is 
positioned to a height just above the lead collar. The 
AGATA detector is then inserted into the frame and po- 
sitioned in the centre of the plate. 

The experimental preamplifier pulse shapes observed 
as a function of the interaction position of a beam of j 
rays can be recorded. The signals are digitized by uti- 
lizing either the GRT4 VME modules [26], which sam- 
ple the signal with an 80 MHz frequency over a 14-bit 
dynamic range, or the GRETINA VME digitizer mod- 
ules, which have a 100 MHz sampling frequency and 
the same 14-bit dynamic range. Both digitizer systems 
utilize an external trigger provided by the core-energy 
signal from the AGATA detector. Data is read out and 
recorded simultaneously for all 37 channels from a sin- 
gle crystal. 

Scan data can be collected in singles and coincidence 
modes. Singles data yields x-y information on the in- 
teraction position of a y ray. Such a measurement is 
relatively quick to perform; however, the z-position in- 
formation has large uncertainties, as it is only defined by 
the segment size. Coincidence scanning utilizes Comp- 
ton scattering to define a single interaction position in 
x-y-z. Such a methodology demands that the y ray 
Compton scatters to an angle of about 90° in the germa- 
nium detector, depositing the remaining energy in co- 
incidence in a collimated ring of scintillation detectors. 
This method is very precise but can result in a very slow 
procedure due to the low coincidence rate between the 
germanium detector and the scintillators. A full char- 
acterisation of a crystal with a grid of 1 mm takes 2-3 
months. A detailed analysis of singles and coincidence 
data for the prototype AGATA symmetric detectors can 
be found in [27, 28]. 

An example of a '^^Cs singles scan of the front face 
of crystal COOl is shown in Fig. 11. Data were recorded 
on a 1 mm grid for 1 min at each position. The system 
was triggered externally with a low-energy threshold of 
~650 keV on the core contact in order to eliminate un- 
wanted Compton scattered or background events. The 
resulting distributions of the intensity of 662 keV full- 
energy depositions confined to a single segment in rings 
1 and 3 are shown in Fig. 1 1 . The x and y axes represent 
the position of the scanning table in a range of ±40 mm. 
The plots clearly show the segmentation pattern of the 
detector and the presence of the coaxial hole in ring 3. 

The distributions of the r30 (10% to 30% of the 
maximum amplitude) and r90 (10 % to 90 %) rise times 
of the core signal are plotted in Fig. 12. In a coax- 



Ring 1 Ring 3 




^0 -20 20 40 -40 -20 20 40 



X (mm) 

Figure 1 1 : (Colour online) Distribution of the intensity of 662 keV 
full-energy depositions confined to a single segment in rings 1 and 3 
of crystal COO 1 . The colour scale gives the number of counts in the 
662 keV peak recorded in 1 min. 



ial n-type crystal, such as those used in AGATA, the 
r30 distribution is dominated by the drift time of the 
electrons towards the core contact, while the 790 distri- 
bution is determined by both the electron and the hole 
transport. As a result, 730 is expected to increase as 
the interaction point is moved from the core contact to 
the outer electrode, while a minimum value for 790 is 
expected at the locations in which the electron and hole 
collection times are equal. This behaviour is confirmed 
by the plots shown in Fig. 12. In ring 3, 730 ranges 
from ~30ns for small radii to ~ 100 ns for large radii, 
while 790 ranges from ~ 100 ns to ~240ns. For ring 1 
(front ring of the crystal) both the 730 and 790 distri- 
butions have a different character due to the more com- 
plex electric field distribution in this region of the de- 
tector. Minimum rise times are observed at small radii 
due to the quasi-planar nature of the electric field lead- 
ing to short charge collection times through the 13 mm 
distance from the front face to the hole drill depth in the 
crystal. These plots also show the influence of the face- 
centred cubic lattice orientation of the Ge crystal on the 
rise times. For pulses measured at the same radius, a 
maximum variation of 30 % with respect to the crys- 
tal axes is observed for the time required to collect the 
charge carriers. This effect must be taken into account 
in the theoretical simulation, if a reUable validation is to 
be achieved. 

The coincidence between an AGATA detector and an 
array of scintillator crystals can be used to select inter- 
actions at a specific location within the crystals. In prac- 
tice, several events for each location are needed in order 
to average the corresponding waveforms and eliminate, 
as much as possible, the effects of the noise. The averag- 
ing procedure, performed for each location, starts with 
a scaling of each waveform by pre-calculated gain fac- 
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Figure 12: (Colour online) [Left] The distribution of the T30 (10- 
30%) rise time of the core signal for 662 keV full-energy depositions 
confined to a single segment in ring 1 (top) and 3 (bottom) of crystal 
COOl. [Right] The corresponding T9Q (10-90%) distribution. The 
crystallographic axes are illustrated. The colour scale gives the rise 
time in ns. 



tors, derived from the Eu baseline difference energy 
calibration, and with a baseline subtraction. The base- 
line is derived for each individual trace from an average 
of the initial 10 samples in each trace. 

Waveforms are then interpolated to allow for more 
accurate time alignment. The pulse amplitudes are sub- 
sequently normalized in order to have the same max- 
imum amplitude for all of them. Finally, the best fit 
average waveforms corresponding to each location are 
obtained through a minimisation procedure. Only 
the central contact, the segment with net charge deposi- 
tion and its neighbours are considered in the fit. "Noisy" 
events which give a large contribution are excluded 
from the fit procedure. The final result is exemplified in 
Fig. 13, where the average (thick red line) and the con- 
stituent (thin blue lines) pulse shapes are shown for a net 
charge deposition in segment c3. The effect of the can- 
cellation of the random noise across the pulses is clearly 
visible. The standard deviation of the baseline noise for 
the average pulses is 0.9 keV, as opposed to 4.7 keV for 
single pulses. The pulse shapes illustrated in Fig. 13 
also demonstrate the signal induced on the neighbour- 
ing segments b3, c2,c4 and d3. These transient signals 
are those induced on adjacent electrodes to the primary 
interaction due to the drift of the charge carriers inside 
the germanium crystal. 
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Figure 13: (Colour online) Average (thick red line) and all constituent 
(thin blue lines) pulse shapes for a typical interaction in segment c3 
following the;j'^ rejection (see text). The signal induced in the core 
and in the neighbouring segments b3, c2, c4 and d3 are also shown. 



4.2. Orsay scanning system 

The Orsay scanning system is based on the same con- 
cept as the Liverpool system, i.e. a well collimated 
strong radioactive source, an accurate moving system 
and an array of scintillator detectors to define the z po- 
sition of the scattering through coincidence measure- 
ments. The main difference is that the z coordinate can 
be continuously scanned due to the absence of the scat- 
ter collimators, which are replaced by special tungsten 
collimators described below. 

A schematic view of the scanning setup based at 
CSNSM Orsay is shown in Fig. 14 [29]. A 477 MBq 
'^^Cs source is encapsulated in a stainless steel cylin- 
drical container with a diameter of 4 mm, a height of 
6 mm and window thickness of 0.4 mm. This container 
is inserted into a collimator made of densimet (W-Ni- 
Fe alloy, density 18.5 g/cm^*). The y rays emitted by 
the source are collimated by a hole with a diameter of 
1.6 mm and a length of 155 mm. 

Six modules of the TOHR (TOmographe Haute Re- 
solution) detector are used to perform the coincidence 
measurements. Each module is made of a stack of 80 
tungsten plates with a triangular shape and a thickness 
of 200 |J,m. These plates have 400 [xm diameter holes po- 
sitioned on a hexagonal lattice and the geometry of each 
stack acts as a many slit collimator (about 8000 slits) 
with a focal distance of ~7 cm. At the back of each stack 
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of plates, there is a Nal(Tl) crystal for the detection of 
y rays. The six modules are positioned in a compact 
semi-circle at +10.2° from the horizontal plane around 
the AGATA detector, all having the same focal point in 
the germanium crystal. The position of the common fo- 
cal point can be changed by translating the TOHR array 
or by turning the AGATA detector about its central axis; 
thus allowing for a full 3D scan of the AGATA detector 
A more detailed description and measurements can be 
found in Ref. [30]. 




AGATA Ge 



bton Scatter 




Figure 14: A schematic diagram sliowing the CSNSM Orsay scanning 
table assembly (top) and a closer view (bottom) showing the 6 Nal(Tl) 
detectors and the mechanical support for the TOHR, the AGATA de- 
tector and the collimated intense '^'Cs source. 



43. Other scanning systems 

Three additional scanning systems are under con- 
struction and validation within the AGATA collabora- 
tion. At the University of Salamanca a system, which is 
similar to the Liverpool and Orsay systems, is being de- 
veloped, while different concepts of scanning are under 
validation in Strasbourg and at GSl. The technique used 
in Strasbourg is based on the pulse shape comparison 
scan (PSCS) principle [31]. In this setup the detector is 



rotated instead of the source in order to scan the front 
and lateral sides of the entire germanium crystal with a 
collimated source of '^^Cs. 

A new characterisation method has been developed at 
GSl [32, 33]. This approach combines the PSCS prin- 
ciple with an advanced y-ray imaging technique based 
on positron annihilation Compton scattering [34, 35]. 
The imaging is accomplished by means of a position- 
sensitive scintillation detector (PSD) [36] coupled to a 
•^■^Na source with an activity of 300 kBq. In order to re- 
duce the background due to scattered y rays, the ^^Na 
source is placed in the centre of a tungsten collimator 
with two conical openings, one of which is oriented 
towards the germanium crystal, the other towards the 
PSD. The setup with the PSD, the source and the col- 
limator can be rotated around a vertical axis at which 
the germanium crystal is located. Two measurements 
are typically carried out, in order to shine and image 
the entire crystal from two different sides. The PSD is 
operated in coincidence with the AGATA crystal. For 
each AGATA-PSD coincidence event, the traces from 
the core and all 36 segments are digitised and stored. 
The position of the y-ray interaction in the PSD is also 
recorded. The latter allows for the 3D reconstruction of 
the trajectories of the two 511 keV annihilation quanta, 
which are assumed to be collinear. By applying the 
PSCS principle, the detector response for a particular 
y-ray interaction point x-y-z in the AGATA crystal can 
be determined by comparison of the two data sets of 
pulse shapes, corresponding to the two scanned sides. 
The duration of a full scan of an AGATA detector on a 
2 mm pitch grid will be strongly reduced compared to 
existing systems, down to about 10 days at Strasbourg 
and 3 days at GSL 

5. The AGATA infrastructure 

The AGATA infrastructure includes the mechanical 
structure and all services to the detectors, mechanics, 
electronics and to the data acquisition system to ensure 
that spectrometer operates reliably. It includes the "life- 
support system" for the detectors providing the cryo- 
genic cooling, the low- and high-voltage power sup- 
plies, constant monitoring, a user-friendly interface and 
reports on critical situations (detector warm up, power 
losses, etc.). This system is called the detector-support 
system (DSS). 

5.7. The AGATA mechanics 

AGATA requires a mechanical structure to accurately 
support the detector elements and enable their safe in- 
sertion and removal. The structure needs to be able to 
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locate the detectors accurately with minimal space be- 
tween each ATC detector in order to maximise the solid 
angle coverage. 

At LNL, AGATA is located at the target position of 
the PRISMA [37] magnetic spectrometer. The main 
design constraint imposed by PRISMA is that both 
AGATA and PRISMA must be free to rotate around the 
beam direction such that the optical axis of PRISMA 
ranges from 0° to 1 17°. Both AGATA and PRISMA are 
therefore mounted on a rotating platform that also sup- 
ports the front-end electronics (digitisers), power sup- 
plies and autofill system. In addition, the arrangement 
at LNL allows the coupling to a range of ancillary de- 
tectors for specific measurements. This setup, including 
all the detector systems, is described in detail in Ref 
[38]. 

The generic support structure for AGATA consists of 
a number of identical flanges (Fig. 15), one for each 
detector module. These are assembled together to pro- 
duce a solid structure as shown in Fig. 16. This generic 
support structure is modular in concept and can be ex- 
panded up to the full 4jt system with 60 such flanges 
and will be used at all the host laboratories. The frame 
to support this structure will be different at each site be- 
cause of the details of the location and the coupling to 
different spectrometers. 




Figure 15: (Colour online) A computer aided design drawing of one 
of the main AGATA support flanges. 

Accurate positioning of the detectors on the support 
structure is achieved by using a system of 3 precision 
sliding rods, one of which is threaded to enabled con- 
trolled insertion and removal. The accuracy of this 
system was checked using a high-precision coordinate 
measuring machine. Under operational conditions the 
endcap underwent some deflection, which was mod- 
elled by finite element analysis (Fig. 17). An adjustment 
mechanism for the detector modules for each flange was 
therefore incorporated into the system. This adjustment 
mechanism comprises 3 rings to provide for the full 6 
degrees of freedom adjustment. In addition, a stand- 




Figure 16: (Colour online) A photograph of 15 flanges mounted at 
LNL. 

alone gauge was manufactured such that the detector 
module could be adjusted to its final orientation before 
mounting into the main support structure. The detector 
module needs to be positioned within 0.1 mm of its the- 
oretical position, and so it was critical that this gauge 
was made to a high accuracy. 




Figure 17: (Colour online) Right: measured endcap deflection of the 
detector module. Left: Finite element analysis deflection calculations 
of the endcap. 



5.2. Detector-support system 

The AGATA DSS consists of the low- and high- 
voltage power supplies, the autofill system and an unin- 
terruptable power supply (UPS) system. The slow con- 
trol of the DSS is managed by a system based on a pro- 
grammable logic controller (PLC), which is accessed 
and controlled from a custom made graphical user in- 
terface. 



15 




Figure 18: (Colour online) Computer aided design drawing of an ATC 
detector and the detector mounting arrangement comprising the main 
flange, 3 precision rods and the detector adjustment rings. 

5.2.7. Low-voltage power supply unit 

This unit suppUes low voltage power to the preampli- 
fiers, high-voltage modules, liquid nitrogen level mea- 
surement system, digitisers and to the PROFIBUS-DP 
(DecentraUzed Peripherals) field bus, which is used for 
controlling the power supplies. 

The power consumption of the digitisers is much 
larger than of the other units. In order to reduce the 
noise pick-up, the digitisers are located as close to their 
power supplies as was practically feasible, at a distance 
of about 8 m. The preamplifiers, high-voltage modules 
and the liquid nitrogen level measurement systems are 
about 15 m from their related power supplies. 

A low-voltage power supply from the company AXIS 
was chosen. It consists of a 4U crate (see Fig. 19), 
which contains all the needed power supplies: +6 V and 
±12 V for the preamplifiers, +6.5 V for the high-voltage 
modules, ±12V for the liquid nitrogen level measure- 
ment system, +48 V and +5 V for the digitisers and 
+24 V for the PROFIBUS-DP field bus. All power sup- 
plies are linear, including the +48 V, which represents 
the most powerful module with a 1440 W load. 

The design of the low voltage supply for the 111 
preamplifiers uses a floating supply system, ensuring 
that the V reference is at the detector itself and all re- 
turn currents pass through the supply to minimise the 
overall detector noise. The preamplifier load is the high- 
est on the +6 V line and, to a lesser extent on the -6 V 
line. The voltage drop across 15 m cable means that 
the +6 V and -6 V voltages have to be regulated at the 
load. This has been achieved to within a strict tolerance 
of ±0.05 V. 




Figure 19: (Colour online) The AXIS low-voltage power-supply unit 
generating ±6 V and ± 1 2 V to power the 111 preamplifiers of the ATC 
detector, +48 V and +5 V for three digitisers, +6.5 V for three HV 
Modules and +24 V for the PROFIBUS-DP network. 

5.2.2. High-voltage module 

Presently, the high voltages for the ATC detectors are 
provided by a standard HV power supply produced by 
the company CAEN. A development is ongoing to equip 
each ATC detector with three HV modules, which will 
generate locally the ~5 kV needed to bias the crystals; 
thus avoiding long HV cables and therefore reducing the 
pick-up noise. In this new design, the high voltages are 
filtered, controlled and monitored directly at the ATC 
detector. The high voltages are produced from a low 
voltage (+6.5 V), which is provided by the AXIS unit. 

Prototypes of the HV modules have been produced 
and they rely on a compact HV unit made by the com- 
pany ISEG. The HV unit is able to communicate with 
the PLC via the PROFIBUS-DP field bus by using cus- 
tomised control electronics. Additionally, a microcon- 
troller provides an autonomous bias shut down opera- 
tion: each HV unit is able to initiate a bias shut down if 
the temperature of the detector or the current measured 
is too high. The bias shut down due to too high tem- 
perature is implemented by assigning one of the three 
HV modules of the ATC detector to act as a master and 
to read the PTIOO temperature gauges of the detector. 
For security and maintenance reasons, several param- 
eters are embedded in the microcontroller: the ID of 
the HV module, the serial number, the self-calibration 
values of the ISEG DC/DC converter, the maximum 
voltage, the maximum ramp-up and ramp-down voltage 
rates and the set threshold value for the current limit. 

A photo of the prototype HV module is shown in Fig. 
20. 

5.2.3. Autofill system 

The liquid nitrogen cooling of the ATC detectors is 
managed and monitored by the AGATA autofill system. 
This system provides information regarding the detec- 
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Figure 20: (Colour online) The AGATA high-voltage module com- 
posed of a compact ISEG DC/DC module and control-electronics 
boards including a microcontroller. 



• A Valve control crate, which contains 4 valve con- 
trol cards, one master card and one dry contacts 
board. The valve control crate is presently under 
development. Parts of its functionality are mim- 
icked using a system based on relays in the present 
setup. 

5.2.4. Uninterruptable power supply 

Many of the AGATA devices are protected by a UPS 
system. The power consumption of the digitisers is 
much larger than of the preamplifiers, which in turn is 
larger than the power consumption of the autofiU sys- 
tem. In case of a power failure, the DSS determines 
when the power supply of each component is to be shut 
down and it initiates the power-off procedure. The goal 
of the DSS is to save enough power to enable a complete 
liquid nitrogen filling cycle. 



tor temperature and the amount of liquid nitrogen in the 
detector Dewar [18]. It controls the filling cycle and op- 
erates all associated valves. 

The autofill system is based on a PLC driven pro- 
cess running under the PROFIBUS data acquisition and 
command standard. It performs the regular filling of the 
detectors, forces the filling of any detector if its temper- 
ature exceeds given threshold, issues warnings if any pa- 
rameter declines from its regular values and issues alerts 
if any of the parameters is beyond the range of noncriti- 
cal values. The system is able to manage a direct filling 
from a pipeline as well as via buffer tanks with the as- 
sociated parameters, warnings and alarms. It can run in 
fully automatic mode, in semi-automatic mode (manual 
start of the fill but automatic follow up of the complete 
fill procedure; remote action using the graphical user 
interface of the GUI), in manual mode (full manual fill 
processed remotely using the graphical user interface of 
the DSS) and in local mode (manual fill of the detectors 
using keys on the autofill hardware in the experimental 
hall). 

Technically the autofill hardware is made of the fol- 
lowing components: 

• A PLC which controls the autofill routine. 

• A PROFIBUS crate, which contains the 
PROFIBUS terminals with various functions: 
PTIOO readout, analogue readout (4 mA to 
20 mA), digital input and output terminals, 
PROFIBUS watchdog, and the relevant power 
supplies. The signals from and the commands 
to the executors (valves, dialers) are sent via the 
valve control crate. 



5.2.5. DSS slow control architecture 

As reliability and safety of the DSS are of paramount 
concern, a PLC has been chosen to run the processes 
and to be the interface to the graphical user interface 
of the DSS. The PLC communicates with the different 
elements of the DSS (autofill, low-voltage power supply 
and high-voltage modules) via the PROFIBUS-DP field 
bus (Fig. 21). 
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Figure 21: (Colour online) PLC architecture of the DSS. The 
PROFIBUS-DP field bus is indicated with red lines. 
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The communication between the PLC and the com- 
puter running the graphical user interface is performed 
through the use of an OPC (object linking and embed- 
ding for process control) server. 

The maintenance server, which is based on an indus- 
trial PC, can be accessed remotely. It gathers all data 
concerning the particular process in order to be able to 
diagnose problems coming from the process itself, from 
the PLC or from the field bus. 

5.2.6. Graphical user interface 

The graphical user interface (GUI) of the DSS is de- 
signed to control and monitor the low voltages, high 
voltages and the autofill parts of the DSS. To ensure 
a user-friendly and easy-to-use interface, different re- 
quirements have been taken into account. Firstly, there 
are different groups of users that demand different kinds 
of information to be displayed and functions to be of- 
fered. Secondly, the amount of data to be displayed is 
quite large. Therefore, it is distributed among different 
display panels, such that only the relevant part is shown 
(Fig. 22). The DSS GUI offers various ways to access 
the whole set of data. For example, there are three tables 
showing all parameters for one detector, user selected 
parameters for all detectors or the raw data blocks from 
the PLC. A graphical representation of historical data 
can be displayed. The user can open multiple tabs and 
for each of them select the items to be displayed and 
updated periodically. 
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Figure 22: (Colour online) The graphical user interface of the AGATA 
detector- support system showing the autofill page. The autofill man- 
ages detectors in groups of 8 units. On the left hand side of the figure, 
the level of liquid in the main supply tank is shown as well as its inter- 
nal pressure and the time remaining to the next programmed fill. On 
the right hand side, the colour of the detector endcap gives an indi- 
cation of the detector temperature (blue for cold, orange for warming 
up, red for warm). Trending (real-time plot of any variable parame- 
ter). Control (parameter setting) and Alarms (alarm status) pages are 
also accessible from this GUI. 



The object-oriented approach allowed the separation 
of the two main threads, the GUI and the PLC commu- 
nication part. Both threads communicate with the DSS 
hardware interface, which holds all the data and man- 
ages the reading and writing from both sides. 



5.3. Cable management 

The management of the many cables required to ser- 
vice, control and extract signal information from the de- 
tectors is very important and has to be taken into ac- 
count in the overall design. The type, length, weight 
and routing of cables has to be considered to enable de- 
tectors to be installed and extracted easily. The cabling 
management system will also be different at each host 
site. 



5.4. Grounding and electromagnetic compatibility 

To provide the best signal-to-noise ratio for the detec- 
tor signals, especially for their use in pulse-shape anal- 
ysis, grounding is of utmost importance. The design of 
the grounding includes the mechanics, the detectors and 
their preamplifiers, the digitisers and all hardware of the 
DSS. The rest of the AGATA system (pre-processing 
electronics and data acquisition system) is optically iso- 
lated from these front-end components. 

It is mandatory that the AGATA grounding system 
shields the various components from both low- and 
high-frequency perturbations. Low-frequency distur- 
bances (~50Hz) cause energy resolution degradation 
in the detectors, but they are generally easy to filter 
out. High-frequency perturbations (from ~lMHz up 
to ~ 100 MHz) are generally more difficult to be fil- 
tered out, and they affect the performance by distort- 
ing the pulse shape of the signals, which might have 
a significant impact on the quality of the pulse-shape 
analysis process. A mesh grounding system is used for 
AGATA with the whole system grounded to a common 
voltage, which is provided by a large common conduc- 
tive plate. The grounding of the mechanical structure 
is formed by interconnected conductive components us- 
ing as short and thick grounding shunts as possible. In 
addition, the power distribution to the front-end com- 
ponents of the array is ensured via a single UPS. The 
measured electromagnetic compatibility performance of 
the AGATA front-end electronics is such that the 50 Hz 
noise is less than 100 (xV RMS and the high-frequency 
noise in the range from 0.1 MHz to 100 MHz is less than 
5mVRMS. 
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Table 3: Count-rate specifications and current limits in kHz after the 
front-end electronics, global trigger and pulse-shape analysis farm. 
The rates are given per crystal and for the setups with 15 and 180 crys- 
tals (AGATA 4;r) with the nominal source-detector distance (23.5 cm). 



System 
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9000 
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3000 (M^ = l)"* 


shape 






300 (My = 30)'* 


analysis 









" With the FWHM value degraded by about 50 % com- 
pared to the nominal FWHM value obtained with a sin- 
gle crystal counting rate in the range 10 kHz to 20 kHz. 

* Limit defined by the check idle cycles. 

'= Current limits of PSA processing. Writing the signal 
traces to disk results in a reduction of the rate to 1 kHz 
to 3.5 kHz per crystal. 

Full AGATA specifications rates with on-line PSA and 
tracking. 



6. Front-end electronics 

The objective of the AGATA front-end electronics 
(FEE) is to digitise the signals from each crystal, pro- 
cess them in real time to establish when the crystals 
detected a y ray (marking the time with a timestamp), 
determine the amount of energy deposited by the in- 
teraction of the 7 ray in each segment and extract the 
total y-ray energy from the core contact of the crystal. 
The positions of the y-ray interactions are calculated in 
the pulse-shape analysis (PSA) farm (see section 10). 
Therefore, the AGATA FEE supplies a short trace of the 
digitised leading edge of each pulse along with the en- 
ergy and the timestamp. The design of the AGATA FEE 
has been a challenge due to the high acquisition rates 
specified for the system. The rate specifications and the 
measured values for the dififerent subsystems are sum- 
marised in Table 3. 

The PSA requires each digitizing sampling ADC to 
be aligned with all the others and that time-stamped data 
based on a common clock are generated. AGATA clock 
distribution and time stamping is linked to a global- 
trigger mechanism, which allows either free running 
time-stamped operation or hardware-triggered opera- 
tion depending on whether data-rate reduction is needed 
in order not to saturate the PSA processor farm. The 
data sent to the PSA contain the same parameters, in- 
cluding timestamps, in both cases. The only difference 



is whether part of the data are rejected prior to the PSA, 
for example because of a multiplicity condition or a re- 
quirement for coincidences with one or more ancillary 
detectors. 
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Figure 23: Schematic view of the AGATA front-end electronics and 
data readout system. 

Fig. 23 shows the schematic design of the AGATA 
FEE and data-readout system, which treats each crys- 
tal (one core plus 36 segment outputs) as a separate en- 
tity. Within that entity the core is treated diff'erently to 
the segments. The core signal is formed by the sum of 
the charge released in all the interactions in the crystal. 
Therefore, it can be used as a trigger for the whole crys- 
tal. The AGATA electronics for each crystal consists of 
the following components: 

• One digitiser comprising 6 segment cards, 1 core 
card, 2 power-supply cards and 2 control cards. 

• Two pre-processing caiTier cards in the Ad- 
vanced Telecommunications Computing Architec- 
ture (ATCA) card format, each containing 4 Com- 
mon Mezzanine Cards (CMC) with PCI Express 
readout to the PSA farm. Seven CMC mezzanines 
correspond to the 6 segment and 1 core card in 
the digitiser and one contains the interface to the 
global trigger and clock system. 

The various elements of the FEE are connected through 
optical fibres in order to achieve appropriate data- 
transmission rates and to maintain good electrical iso- 
lation. 

6.1. The digitiser 

The principal goal of the digitiser is to interface the 
detectors with the AGATA signal-processing system. In 
order to do that, the digitiser module performs the fol- 
lowing tasks: 
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• Receives all 37 preamplifier outputs from one crys- 
tal. 

• Digitises the input signals at a rate of 100 MHz, 
using 14-bit ADCs. 

• Serialises the ADC data and transmits it over opti- 
cal fibres to the pre-processing electronics. 

• Implements a constant-fraction discriminator 
(CFD) algorithm in an FPGA to generate an 
isolated fast-logic signal to be used in hardware 
triggers of ancillary detectors. 

• Implements a time-over-threshold (TOT) algo- 
rithm using preamphfier inhibit signals (see section 
3.3). 

• Provides spare digitiser channels and inspection 
lines for maintainability and diagnostics. 

• Provides interfaces for re -programming, control in 
an electronics and detector workshop, as well as an 
interface for slow control. 

The AGATA specifications require that the digitiser 
module is mounted less than 10 m from the detector's 
preamplifiers to minimise signal degradation. The digi- 
tiser is housed in a water cooled box of the size 30 cm 
X 14 cm X 55 cm and it contains two modules with the 
electronics for two crystals. A block diagram of the 
digitiser is shown in Fig. 24 and a photograph of a 
printed circuit board for one module is shown in Fig. 25. 
The main building blocks of the digitiser are: the differ- 
ential analogue input, the flash analogue to digital con- 
verter (FADC) sampling blocks, the field programmable 
gate arrays (FPGAs), which receive the FADC data and 
serialise it, the optical transceiver blocks, the global 
clock receiver block and the monitoring block. The 
monitoring block provides control of the spare channels 
and inspection lines via a slow control link to the GUI 
based system control software. The temperature of the 
various parts of the digitiser electronics is also acces- 
sible by the slow control and a local over-temperature 
shutdown is implemented in the control/power supply 
card. 

The analogue input buffer adds the optional off'set sig- 
nal and also includes an anti-aliasing filter before the 
FADC. The FADCs convert the signals into 14 bits at 
100 MHz; 15 data bits ai-e sent to the FPGA for serial- 
isation and transmission as 16-bit data over the optical 
fibres. Bits to 13 are FADC data, bit 14 is the FADC 
overrange and bit 15 is used as the synchronisation pulse 
from the GTS (for a description of the GTS, see subsec- 
tion 6.3). 
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Figure 24; Block diagram of tlie AGATA digitiser. 



Although the main purpose of the FPGAs in the seg- 
ment and core cards is to serialise data, the FPGA of the 
core card also implements two algorithms for real-time 
data processing. The first is a digital CFD, which is used 
as a trigger for ancillary detectors (output both on an 
isolated connector and over a fibre). The CFD algorithm 
runs in the digitiser rather than in the pre-processing be- 
cause of the latency of the fibre links (several 100 ns) 
and the need for a prompt output. The second algo- 
rithm is a TOT algorithm that measures the time during 
which the preamplifier inhibit signal is asserted when it 
receives an over-range energy input. Since the pream- 
plifier recovers from overload by discharging a capaci- 
tor using a constant current source the period of over- 
load (TOT) is a measure of the energy in the crystal 
when the preamplifier is over-range (see section 3.3). 
A special serial fibre protocol is used to indicate when 
FADC data is replaced by TOT data. 

One 8 channel fibre-optic transceiver link (4 RX, 4 
TX) is used to interconnect the core FADCs serialised 
data and the timing and control of the pre-processing 
subsystem. The "transmit" channels are allocated to the 
core 100 MHz FADC data (two channels), the core CFD 
logic signal (one channel), and the global clock calibra- 
tion feedback path (one channel). The "receive" chan- 
nels are used for the global clock signal, the synchro- 
nisation signal, and the analogue offset control for the 
core electronics. A 12 channel multi-fibre transmits the 
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Figure 25: (Colour online) Photograph of the AGATA digitiser. 



6 serialised segment FADC data streams from each of 
the segment cards. Each analogue channel from the de- 
tector can be inspected before and after the digitisation. 

One analogue line is provided before and two after 
coding per group of 6 segments channels. The global 
clock function reconstructs a high quality clock with a 
very small jitter (<7 ps) from the clock received through 
one optical receiver within the core transceiver 

The digitiser is able to work in stand-alone mode 
when a specific firmware and GUI interface (called 
SAMWIZE) are used. In this mode it is possible to vi- 
sualise the recorded traces from the detector input sig- 
nal (oscilloscope mode). It is also possible to record up 
to 6 traces per segment and one per core. The Multi- 
Channel Analyzer mode is capable of histogramming 
up to 6 segment channels and one core in parallel. This 
mode was particularly useful for developing and test- 
ing the core energy and constant-fraction discriminator 
(trigger) block of the pre-processing firmware. 

6.2. The pre-processing electronics 

The pre-processing system reduces the data volume 
from the digitisers by a factor of about 100 by extract- 
ing and processing data from the digitiser's data stream 
only for the segments which have registered a detection 
of a y-ray interaction. Further filtering (triggering) can 
optionally be performed in conjunction with the GTS 
from which a clock for the digitiser and the timestamp 
information is also derived. The pre-processing sends 



the filtered data to the PSA farm. The processing rate 
for traces in the segments is the same as in the core. 
The triggering is always made by the core, so the core 
contact electronics is the master and the segments are 
controlled by it. 

The pre-processing hardware de-serialises and pro- 
cesses the incoming data streams and stores the traces. 
The sampling speed is 100 MHz in the reconstructed 
data streams, so the pre-processing hardware also uses 
100 MHz clock rates for incoming data (some internal 
clocks are running at 200 MHz). The GTS interface pro- 
vides the system clock and a trigger system. 




Figure 26: Block diagram of the AGATA pre-processing mezzanines. 




Figure 27: (Colour online) Photograph of the AGATA pre-processing 
ATCA carrier card. 



The AGATA trigger system can be used to reduce the 
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Figure 28: (Colour online) Photograph of the AGATA pre-processing 
mezzanine card. 



counting rate. Where rate reduction (see Table 3) is not 
required, the pre-processing runs in triggerless mode, 
which means that all the processed data are sent to the 
PSA farm. In this case a software trigger is performed 
after PSA and tracking. The maximum delay (latency) 
which can be accommodated in the pre-processing hard- 
ware while waiting for a trigger decision is limited by 
the maximum of the trace length, which is 20 [xs. Lower 
values can be configured in the trigger system, but coin- 
cidences separated by more than 20 (,is are detected by 
using software triggers. 

The local trigger of the core signal is a digital-trigger 
algorithm operating on the data stream from the core 
contact of the crystal. When this trigger finds a pulse 
in the data it generates a local trigger output which in- 
dicates to all the segment electronics that they should 
also extract a trace from the data stream. Traces are 
stored locally, within each pre-processing channel. The 
traces are held in each charmel's local memory for up 
to 20 yLS while the GTS makes a decision. An event can 
be either accepted or rejected. For events which are ac- 
cepted, the pre-processing stores a trace of the digitised 
leading edge of the pulses from the core and from all 36 
segments in a buffer waiting to be sent to the PSA farm. 

In addition to selecting useful portions of the in- 
coming data stream using a trigger algorithm, the pre- 
processing also applies the moving window deconvolu- 
tion (MWD) algorithm [39, 40] on the incoming data 
streams to determine the y-ray energy deposited in each 
segment and in the whole crystal. The MWD filter out- 
put is sampled a programmable time after the trigger 
algorithm detects the start of a digitised pulse in the in- 
coming core data stream. The energy parameter sam- 
pled from the MWD filter is stored along with a data 
trace showing the samples taken during the leading edge 



of the pulse and a timestamp indicating when the trigger 
happened. 

Core and segment mezzanines have PowerPC proces- 
sors with an embedded Linux operating system hous- 
ing all the drivers and applications for the card con- 
trol and online temperature monitoring. A slow con- 
trol GUI, based on a web server architecture, has been 
developed for the control and online monitoring of the 
status of the mezzanine cards. The block diagram of 
the pre-processing mezzanine is shown in Fig. 26. Pho- 
tographs of the printed circuit boards of the carrier and 
mezzanine cards are shown in Fig. 27 and 28, respec- 
tively. Carrier cards, implemented in the ATCA stan- 
dard [41], house the CMC format mezzanines, which 
read data from the core and segment cards in the digi- 
tiser and a GTS mezzanine, which communicates with 
the GTS system. 

A block diagram of the GTS mezzanine is shown in 
Fig. 29 and a photograph of the printed circuit board in 
Fig. 30. The GTS mezzanine receives the global clock, 
aligns it locally within the pre-processing and then uses 
a dedicated data path in the core mezzanine to align the 
clock in the digitiser It is also used to accept or reject 
local triggers. The core and segment mezzanines re- 
ceive data from the digitiser and process the data using 
digital filter algorithms. The mezzanines are read out 
via the carrier, all data being concentrated into a single 
FPGA per carrier before transmission via PCI express 
to the PSA farm on demand by the PSA. The carrier 
reads each mezzanine at lOOMB/s. A different FPGA 
handles the trigger interconnections within the carrier. 
Trigger connections between the two carriers handling 
one crystal use a dedicated backplane link (TCLK). The 
ATCA carrier cards accept 4 CMC format mezzanines. 
The connection between the CMC and the carrier card is 
achieved by two Mictor connectors, each with 1 14 pins. 

6.3. The Global trigger and synchronisation ( GTS) sys- 
tem 

Data synchronisation is an important aspect in the op- 
eration of the trigger and readout systems of AGATA. 
Tracking and PSA require the concurrent digitisation 
of preamplifier signals of the 36 segmented Ge crys- 
tals composing the array. Therefore, the design of the 
front-end readout and level- 1 (LI) trigger in AGATA 
follows a synchronous pipeline model: the detector data 
are stored in pipeline buffers at the global AGATA fre- 
quency, awaiting the global LI decision. The LI latency 
must be constant and match the pipeline buffer length. 
The whole system behaves synchronously and synchro- 
nisation at different levels and in different contexts has 
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Figure 30: (Colour online) Photographs of the two sides of the GTS 
Figure 29: Block diagram of the GTS mezzanine. mezzanine. 



to be achieved and monitored for proper operation of the 
system. 

In AGATA each crystal is considered as a separate 
entity and, from the point of view of the data acquisi- 
tion system, the whole detector may be considered as 
the aggregation of synchronised data supplied by indi- 
vidual crystals, possibly disciplined by a global trigger 
primitive. 

The data from the core contact are processed for event 
detection and either an LI trigger request or a local 
trigger is generated. The choice between the two be- 
haviours is done upon configuration, the former corre- 
sponding to an effective way to reduce front-end data 
rates in cases where any one of the stages of the readout 
chain is unable to perform at the actual data throughput. 
From the logical description of the front-end operation 
given above it turns out that a certain number of global 
time-referenced signals are needed. The most impor- 
tant of these are: the common clock, the global clock 
counter, the global event counter, the trigger controls, 
the trigger requests and the error reports. In AGATA, 
the transport medium of all these signals is shared by 
use of serial optical bi-directional links connecting the 
FEE of each crystal with a central global trigger and 
synchronisation control unit in a tree-like structure; thus 
actually merging together the three basic functionalities 
of synchronisation distribution, global control and trig- 
ger processing. 

The common clock is a 100 MHz digital clock sup- 



plied by a central timing unit and used to clock the 
high-speed optical transceivers reaching the FEE of ev- 
ery crystal. At the crystal receiving side the clock is 
reconstructed and filtered for jitter. The clock signals 
of each crystal may be equalised for delay and phase; 
thus accounting for different fibre lengths and different 
crystal locations in the array. 

The global clock counter is a 48-bit digital pattern 
used to tag event fragments before front-end buffer for- 
matting. The pattern is the actual count of the global 
clock. It will be used by the PSA processing and by the 
global event builders to merge the event fragments into 
one single event. The global event counter is a 24-bit 
digital pattern used to tag event fragments before front- 
end buffer formatting. The pattern is the actual count of 
the LI validations. The trigger control must guarantee 
that subsystems are ready to receive every LI accept de- 
livered. This is essential in order to prevent buffer over- 
flows and/or missed trigger signals when the crystals 
are not ready to receive them. In either case, the con- 
sequence would be a loss of synchronisation between 
event fragments. 

The trigger requests are generated by the core contact 
signal from the AGATA detectors by asserting a signal 
that is transmitted via the high-speed serial links of the 
GTS network upwards to the central trigger unit. All 
the trigger requests collected from the crystals at each 
global clock cycle form a pattern that can be processed 
centrally for multiplicity or coincidence with ancillary 
detectors. The result of this processing stage consti- 
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tutes the LI validation. Error reports indicate abnormal 
conditions such as buffer overflows, local faults, built-in 
self tests, etc., and can be reported centrally for proper 
corrective actions. 



GTS global clock 




XTXjXjXIX 



Figure 3 1 : (Colour online) Illustration of the topology of the GTS 
tree. Any detector node can produce at any time a trigger request, 
which is labeled with the clock and transmitted upwards in the tree, 
through the fan-in-fan-out architecture to the global trigger processor. 



The AGATA GTS has a tree topology as shown in 
Fig. 31. The tree originates from the root node, which 
at the same time acts as the source of all global informa- 
tion (clock, timestamps, commands, LI validations) and 
all the trigger requests. It also performs fast monitoring 
of signals and services requests coming from the crys- 
tals. To solve the problems of building a bi-directional, 
high-capacity and high-speed network tree that drives 
hundreds of nodes which are displaced by several tens 
of meters, a certain number of technological issues have 
been addressed. Among these issues are the fan-out of a 
source synchronous transmission, noise immunity, low 
error rate and throughput. The GTS tree is composed 
of five different elements: the root node, the backplane, 
the fan-in-fan-out nodes, the fibre connections and the 
mezzanine interface. 

At the root node, a fast optical link forwards all the 
collected trigger requests to the trigger processor, which 
is a high-performance and fully-pipelined custom-built 
processor that sorts incoming requests, computes mul- 
tiplicities and prompt or delayed coincidences among 
seven user defined trigger partitions. The trigger pro- 
cessor is fully programmable in the C programming lan- 
guage by the user through an application programming 
interface and it is controllable via a standard PCI ex- 
press interface. 



7. Coupling of complementary instrumentation 

The coupling of large y-ray detector arrays, such as 
EUROBALL [4] and GAMMASPHERE [5], to com- 
plementary detectors has played a major role in spec- 
troscopic investigations in the last decade. The AGATA 
array will be used with radioactive or high-intensity sta- 
ble beams and in most experimental conditions, in order 
to fully exploit its capabilities, it will be essential to cou- 
ple it to devices providing complementary information. 
Additionally, the foreseen use of AGATA at different fa- 
cilities makes some of the complementary devices (e.g. 
beam-tracking devices in the case of fragmentation fa- 
cilities) absolutely necessary for the normal operation 
of the array. 

Within the AGATA project, dedicated new comple- 
mentary detection systems will be developed in addi- 
tion to adapting existing ones. The use of fully digital 
sampling electronics is desirable in the future. How- 
ever, in order to utilise the many existing systems based 
on conventional analogue readout electronics an inter- 
face to the GTS system has been developed. This inter- 
face is called AGAVA (AGATA VME Adapter). Since 
nearly all of the presently available detectors within 
the AGATA community have front-end/readout systems 
based on the VME or VXI standard, the AGAVA inter- 
face has been developed in the VME standard with full 
compatibility with the VXI readout modes. 

7.1. General description of AGAVA 

The AGAVA Interface is a one-unit wide A32/D32 
type VME/VXI slave module. It is the carrier board for 
the GTS mezzanine card used in AGATA for the global 
trigger and timestamp distribution. 

The main purpose of the AGAVA interface is to merge 
the AGATA timestamp-based system with conventional 
readout based on VME or VXI, which is used for exam- 
ple by EXOGAM [42]. AGAVA also has the necessary 
connectors to interface with the VME Metronome and 
SHARC [43] link systems. The logic process is con- 
trolled by an FPGA of the type Virtex II Pro. The block 
diagram of the AGAVA interface is shown in Fig. 32. 
The AGAVA module (see Fig. 33) includes all neces- 
sary connections for the trigger cycle and for a total data 
readout system [43]. It contains also a passive Ethernet 
interface, which provides a direct connection to the GTS 
mezzanine card. 

7.2. The operation of AGAVA 

The AGAVA module supports VME access, provides 
interfacing with the GTS mezzanine and provides in- 
put/output signals at the front-panel Lemo connectors. 
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Figure 32: Block diagram of tlie AGAVA interface. 



Two VME access modes have been implemented in 
AGAVA: the VME single-access read or write opera- 
tions, according to the standard VME handshake rules, 
and the chained block transfer (CELT) readout mode 
on the VME bus. Standard access can be used both 
for board configuration and readout, while the CELT 
mode (available in many commercial VME ADC and 
TDC front-end modules) has been introduced to speed 
up data readout. The CELT mode makes use of a com- 
mon address specified during the initialisation phase in 
a dedicated register for all participating modules in the 
chain. The VME CPU reads their data in blocks of vari- 
able lengths starting from the module defined as the first 
in the chain, followed by one or more intermediate mod- 
ules and finished by the one configured as the last one. 
The CELT mode is used for faster readout and is helpful 
for event building. The CELT readout mode of AGAVA 
has been extensively tested. The AGAVA interface also 
implements a VXI EXOGAM-like readout. 

The interface to the GTS mezzanine provides the full 
trigger-request/validation/rejection cycle and acquires 
also the clock tag and the event number, the latter only 
in case of a validated cycle. 

The AGAVA module receives, through the front- 
panel input, the trigger requests from the complemen- 
tary detector (named here external trigger) that is passed 
as a single 10 ns wide pulse to the GTS system (an in- 
ternal trigger request can be generated for testing pur- 
poses). After the trigger request is passed to the GTS 
system, the AGAVA module receives the local trigger 
signal and tag and waits for the validation or rejection 
signal and tag coming from the GTS System. The la- 
tency time to receive the validation or rejection depends 
on the GTS tree architecture and is of the order of 12 ^s. 



Figure 33: (Colour online) Photograph of the AGAVA module with 
the GTS mezzanine card. 



Information is stored in registers or RAM memory 
(depending on the firmware functionality mode). Once 
the information is stored, the AGAVA module sets the 
data ready flag to inform the VME or VXI system that 
the data can be readout by the CPU module and trans- 
feiTed to the event builder. The busy flag is set after 
accepting every trigger request. The release of the busy 
(i.e. ending the dead time for new triggers) strongly de- 
pends on the VME/VXI readout mode and speed. In 
case of high trigger rates, the backpressure front-panel 
input can be used to reduce the trigger request to the 
GTS system. 

In order to support the wide range of complementary 
instruments, two different AGAVA functionalities have 
been introduced: common dead-time mode and parallel- 
like mode. 

In the common dead-time mode only one trigger at 
a time is accepted. The trigger request, acknowledged 
only if the AGAVA module is not in busy status from 
the previous cycle, sets the busy state until all data are 
readout by the CPU module. 

In the parallel-like mode, when a trigger request ar- 
rives on the front-panel input, the AGAVA module sets 
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and holds its busy status only until the local trigger and 
local trigger tag are received from the GTS mezzanine 
card and are stored in the AGAVA multi-event RAM. 
A new trigger request can be accepted afterwards by 
the AGAVA module avoiding the GTS latency time and 
with a sizable reduction of the total system dead time. 
The flow of validation or rejection triggers and tags re- 
ceived by the AGAVA module are compared with the 
local ones present in the multi-event RAM and if they 
match correctly the VME/VXI CPU can read the rela- 
tive data. 

8. Data acquisition 

As shown in Fig. 34, the data acquisition (DAQ) sys- 
tem receives the preamplifier traces from the FEE and 
processes them into several stages up to the storage of 
reconstructed events; see section 8.1 for a description 
of the data flow. Services described in section 8.2 are 
needed to control and monitor the whole system includ- 
ing the electronics. For a complex instrument such as 
AGATA, large computing, network and storage capabil- 
ities are needed, as will be discussed in section 8.3. 




Storagi 



Figure 34: Layout of the AGATA data acquisition system. 



8.1. Dataflow 

The data flow integrates the algorithms needed to pro- 
cess the information from the interaction points. At the 
end of the chain, data have to be provided to the users in 
the AGATA Data Flow (ADF) library [44]. After the 
PSA, data from all crystals are merged together tak- 
ing into account the physics correlations provided by 



the pre-processing electronics. These correlations can 
be based on timestamp or event number depending on 
which AGATA trigger mode is used (see section 6) If 
additional detectors are coupled to the AGATA array, 
the corresponding data flow is assembled in the same 
way. After this procedure tracking is performed to re- 
construct the y-ray trajectories. This procedure is then 
followed by the data storage on a large local disk array 
before being sent to Grid Tierl computing centres based 
at INFN-CNAF in Bologna (Italy) and at CC-IN2P3 in 
Lyon (France). It is also possible to write data to disk at 
each step of the data flow. 

Because of the high rate specifications (50 kHz per 
crystal, see Table 3) and the required processing power 
and large data transfer bandwidths, the algorithms have 
to be distributed over several computers (CPUs and 
cores) At this high rate the storage of 200 B per seg- 
ment and core gives a data rate of about 400MB/s per 
crystal at the output of the FEE. The data flow and em- 
bedded algorithms are managed by the Narval [45, 46] 
framework. Narval is a highly-distributed data acquisi- 
tion system running across a network and acting like a 
single program that transports the data to the storage of 
the reconstructed events. Narval is based on actors cor- 
responding to separate processes that receive and send 
out data at any stage of the data-flow chain. It is devel- 
oped in Ada for a high reliability and safety level and 
it is very flexible since the actors can load C and C++ 
shared libraries. 

The event-builder and merge algorithms are written 
in Ada while the PSA and tracking algorithms are writ- 
ten in C++. The originality of the DAQ system re- 
lies on the development of the ADF framework [44], 
which has been designed in order to facilitate the sep- 
aration between data transport and algorithm develop- 
ment. ADF provides the interface between Narval and 
the algorithms (PSA and tracking). It also allows to en- 
code/decode the data at various stages of the data flow. 
ADF can be either used with Narval or as a stand-alone 
library for data replay and it is part of the GammaWare 
package [47] (see section 12). 

In order to survey the data acquisition system, Narval 
provides a "spy" mechanism by which it is possible to 
sample pieces of the data flow from each actor for con- 
trol, monitoring and online analysis purposes. The spy 
is available to any client connected to the DAQ services 
network. 

8.2. Services 

A set of services are provided to control and moni- 
tor the whole system including the detector electronics. 
The run control centralises information from different 
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subsystems of AGATA. It also collects messages, errors, 
and provides an ELOG logbook [48] for the users. A 
GUI is provided for the electronics and run controls. 

8.2.1. Run control 

The main purpose of the run control (RC) is to control 
and monitor the DAQ components. It coordinates the 
large number of activities that are necessary to achieve 
an operational state for the AGATA detector including 
the DAQ system. Actions like initialisation, setup of 
several components, start and stop of the data acquisi- 
tion, are performed by the operator through the RC sys- 
tem. The RC interacts with the electronics control sys- 
tem, which is in charge of the control of the hardware 
devices. The RC also provides the monitoring of the 
data acquisition (input rate, buffer occupancy, error rate, 
etc.), error report, error handling functionalities and log- 
ging capabilities. The RC is based on the middleware 
produced by GridCC [49]. To communicate with all 
the controlled components a standard WSDL (Web Ser- 
vices Description Language) interface has been defined. 

8.2.2. The Cracow graphical user interface 

The GUI known as "Cracow" [50] allows the user to 
perform commands such as loading the configuration, 
starting and stopping the DAQ system, etc. All the pos- 
sible states of the DAQ state machine are reflected in 
the GUI by a set of buttons that allow particular ac- 
tions to be performed. Cracow is an independent tool, 
which communicates with the RC through the web ser- 
vices. It can be run on any computer connected to the 
AGATA network. As mentioned before, Narval is a sys- 
tem whose actors (algorithms) are highly distributed, 
which means that the user can create actors on many 
different processors. A configuration of these actors is 
called a topology. The Cracow GUI can show a work- 
ing diagram of the currently loaded topology. The graph 
is constructed based on information delivered online by 
the actors. They send information about themselves and 
about their direct collaborators. 

8.2.3. Spectrum server and viewer 

The different algorithms implemented in the DAQ 
provide spectra or matrices. In order to collect his- 
tograms from the various sources, the GRU spectrum 
server (GANIL Root Utilities) [51] has been imple- 
mented. As soon as a root histogram is added in the 
GRU database, it is immediately available on the net- 
work and can be displayed by any client like ViGRU 
[51] or the "Cracow" GUI described above. 



8.2.4. Electronics control 

The electronics control (EC) system is an ensemble 
of software tools that allows to setup of all the elec- 
tronics subsystems (digitisers, pre-processing, GTS and 
possibly ancillary electronics) and to monitor some key 
parameters. The main functionalities to be provided by 
the EC system are the following: the description of the 
system insuring its coherency, the localisation of the 
servers controlling each subsystem, the initialisation of 
the different devices with the correct values, the saving 
of all the setup parameters for the whole electronics or 
a part if necessary, the restoration of a previously saved 
setup, the monitoring of key parameters in the different 
boards, the handling of error/alarm events passing them 
to the RC. The EC system has to insure the integrity of 
the system and the synchronisation of the state machine 
of the different subsystems. It accepts a set of simple 
commands from the RC (setup, go, stop, get state, etc.). 

The EC system is designed with a client/server ap- 
proach and with partitions accordingly to each type of 
hardware. There is one EC subsystem with its engineer- 
oriented GUI for each type of hardware, i.e. one EC 
subsystem for digitisers, one for carrier boards, one for 
core/segment mezzanines, one for GTS and one for an- 
cillary electronics. In order to centraUse the different 
electronics subsystems, a Global Electronics Control 
(GEC) based on the ENX [52] framework has been de- 
signed. The GEC acts as a coordinator between the dif- 
ferent subsystems and the RC. 

8.3. Hardware implementation 

The system first implemented for AGATA at LNL 
has been designed to be scalable to cope with the full 
array specifications, easily movable between different 
host sites and easily maintained by the collaboration. 
The AGATA DAQ hardware (servers, disk, local net- 
work, etc.) and software should be seen from the host 
laboratory as a black box. 

On the hardware side, there is a big requirement for 
computing power to run algorithms hosted by the data 
acquisition (PSA, event-builder, merge, tracking). Each 
crystal is attached to one PSA server. Since the PSA 
processing is the most demanding part of the system, 
Narval distributes events on the different cores of each 
PSA server. There is also one server for event building, 
one for merging with ancillary devices, one for track- 
ing, three disk servers, two for slow control, two for 
DAQ box control, and two for data analysis. All servers 
are identical and connected to a KVM (Keyboard Video 
Mouse) switch to ease the hardware maintenance. The 
model chosen is the IBM® x3550 lU server (from ini- 
tial to M3 versions) with 2 quad cores Intel® Xeon® 
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CPUs each, which allows the PSA algorithms to be run 
at a crystal rate of about 4 kHz. The Debian distribution 
[53] of the GNU/Linux operating system is installed on 
all machines. 

AGATA generates a large data flow, which must be 
stored at very high rate, especially when traces are 
stored. The local disk array of AGATA consists of 4 disk 
servers. Three of them are attached to a SUN® Storage 
Tek^'^ 6540 fibre channel system with 112 TB of disk 
(90 TB using RAID 6) and shares data using the GPFS 
(General Parallel File System) clustered file system. A 
fourth server attached to the storage system provides the 
interface to the Grid for further data storage or analysis 
(see section 12). Many services are needed to handle 
the DAQ box: monitoring, documentation (Wiki doc- 
umentation mirrored on the main DAQ web site [54]), 
network control, installation server, hardware monitor- 
ing using Zabbix [55], etc. These services are hosted 
on a virtual Xen [56] machine. A backup system is in- 
stalled on a SUN® Fire^'^ X4560 system with 8 TB of 
disks to save periodically the vital elements of the DAQ 
box. 

In order to reduce the dependency particular to the 
host laboratories, basic network services are provided 
in the DAQ box: DHCP (Dynamic Host Configuration 
Protocol), DNS (Domain Name System), software dis- 
tribution, etc. External users (developers, system ad- 
ministrators) can access the DAQ box through a VPN 
(Virtual Private Network). 

Networking is an important issue having in mind the 
high data rates foreseen in forthcoming AGATA experi- 
mental campaigns. The communication between nodes 
is organised around two different networks, one for data 
acquisition called DAQ network and the other for net- 
work services called DAQ services. 

The data acquisition infrastructure can be hosted in 
one or several buildings. For example at LNL, the Tan- 
dem accelerator building hosts the data flow and DAQ 
services machines while the Tier 2 computer room hosts 
all the disk servers. Both rooms are connected with 
3 GB optical links. 

8.3.1. Algorithm for compression of traces 

In order to reduce the size of the stored traces a sim- 
ple algorithm has been devised to compress the traces 
to about one half of their original size. The compres- 
sion algorithm is based on the observation that the dif- 
ference between the measured voltages of two consec- 
utive sampling points is small. Fig. 35 shows a typical 
spectrum of differences between neighboring sampling 
points. The spectrum is broad with a very sharp peak 
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Figure 35: (Colour online) Typical difference spectrum of recorded 
traces. The parameter on the horizontal axis is the difference in the 
value of two consecutive sampling points. The shaded region, which 
marks difference values from -31 to +32, contains 93 % of the data. 
The histogram contains data from about 13 million sampling points. 

around zero difference. The mean value of the distri- 
bution is positive because the traces are recorded at the 
beginning of the pulse to cover the leading edge of the 
signal. 

The sampling points have a 14-bit resolution and are 
stored by the pre-processing electronics in 16-bit un- 
signed integer values. The trace length is typically 100 
points. In Fig. 35 one may see that about 93 % of the 
difference spectrum is contained in the region between 
-31 and +32. This fact will be used for the compression 
as several of these small differences can be packed into 
16 bits. The chosen storage format is the following: 

• The first sampling point of the trace is always 
stored with its original value. It becomes the first 
reference point. 

• If three consecutive sampling points after the ref- 
erence point differ from their respective previous 
sampling point by values between -15 and +16, 
the three differences are encoded in 3 x 5 bits in 
a 16-bit unsigned integer with the most significant 
bit set to 0. 

• If one sampling point after the reference point dif- 
fers from the reference by a value between -3 1 and 
+32, the difference is encoded in an 8-bit unsigned 
integer with the two most significant bits set to the 
binary value 10. 

• Otherwise, the 14-bit value is stored in a 16-bit un- 
signed integer with the two most significant bits set 
to the binary value 11. 
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• Finally the reference point is advanced to the last 
point stored so far. 

For most of the sampling points of a trace, three con- 
secutive points from originally three 16-bit words can 
be packed into a single 16-bit word. For traces with 
100 points excluding the headers of the data, a com- 
pression of about 38 % is achieved. The compression of 
data containing both sampling points and data headers 
is typically about 44 %. 

9. Pulse-shape simulations 

In order to realise the real-time localisation of the 
scattering sequence following a y-ray interaction inside 
an AGATA detector, the experimentally digitised pulse 
shapes will be compared with a basis data set. For on- 
line pulse shape analysis (PSA) to be implemented suc- 
cessfully, it has been estimated that more than 30000 
basis sites per crystal for a 2 mm Cartesian grid [57] are 
required. The basis can be obtained from either exper- 
imental or simulated data. However, presently the time 
prohibitive nature of the experimental methodology (see 
section 4), means that the only practical way to generate 
a full basis is by calculation. 

The AGATA collaboration has therefore developed 
electric field simulation codes such as Multi Geometry 
Simulation (MGS) package [58], the Java AGATA Sig- 
nal Simulation (JASS) toolkit [59] and the IKP Detec- 
tor Simulation and Optimisation method now called the 
AGATA Data Library (ADL) [60-63], in order to facili- 
tate the realisation of the on-line PSA. These codes have 
been used to generate variants of a single crystal basis 
and an earlier code originally developed for the MARS 
project [11] was used to optimise the physical segmen- 
tation scheme in the depth of the AGATA crystal. 

In general, the electric field simulation codes use 
finite-element methods for solving partial differential 
equations, such as the solution to Poisson's equation. 
Environments such as FEMLAB and DIFFPACK can 
also be adapted for providing an analytical solution to 
the complex electric field distribution inside the AGATA 
crystals. The codes utiUse the approach shown in block 
format in Fig. 36. The calculations are performed on a 
user specified 3D grid that maps a given detector vol- 
ume. Results from each stage of the calculations are 
stored in matrices. The values at each point in the ma- 
trices are then recalled to generate the pulse shape re- 
sponse, as determined by the trajectories of the charge 
carriers through the weighting field [64]. 

The MGS package [58] utilises MatLab's matrix en- 
vironment to derive the expected pulse-shape response 



at the contacts of any geometry of HPGe detector. The 
stand-alone package has been compiled for use with 
both Linux and Microsoft Windows operating systems. 
MGS was developed as an alternative specialised solu- 
tion to the commercial packages. 
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Figure 36: Data flow diagram for the simulation of the expected pulse 
shapes at the contacts of any ai'bitraiy HPGe detector geometry. For a 
given detector, the crystal volume can be divided into a cubic matrix of 
lattice sites. Values for the electric potential, electric field and weight- 
ing field are calculated at each position. The drift velocity matrices 
are calculated from the electric field matrix. The detector response for 
a given interaction site is calculated by tracking the trajectory of the 
charge carriers through the weighting field [64]. 



The top panel of Fig. 37 shows the geometry of the 
AGATA symmetric prototype crystal as generated by 
MGS. The middle panel shows the electric potential and 
the bottom panel the electric field values in the z-x plane 
at y = 51 mm (centre of the crystal). The results show 
the decrease in potential and field strength as a function 
of increasing radius from the central anode. The figure 
clearly illustrates the complex nature of the electric field 
within the closed-end coaxial geometry of the AGATA 
crystal. 

At high electric fields and low temperatures the 
charge carrier drift velocities in germanium become 
anisotropic [65]. They depend on the electric field vec- 
tor with respect to the crystallographic lattice orienta- 
tion. The electron drift velocities, Ve, are saturated at 
field strengths >3000 V cm"' for the <100> and <1 10> 
directions, and at ~4000Vcm ' for the <111> direc- 
tion. The hole drift velocities vj, are saturated for fields 
~4000Vcm"' along all three major crystallographic 
axes. The models for anisotropic drift of the electrons 
and holes in n-type HPGe detectors implemented in 
MGS have been derived from work published in [60] 
and [66], respectively. 

Before the induced current at each electrode can be 
calculated, the weighting potentials and weighting fields 
must be generated. The calculation is performed with a 
null space charge density and with +1 V on the sens- 
ing electrode with all other electrodes grounded. For 
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Figure 37: (Colour online) An MGS simulation of the AGATA sym- 
metric prototype crystal geometry (top). Electric potential (middle) 
and electric field values (bottom) in the z-x plane at y = 5 1 mm (cen- 
tre of the crystal). The results show the decrease in potential and field 
strength as a function of increasing radius from the central anode. 



a given interaction position, the charge-pulse response 
observed at any electrode depends on the trajectory of 
the charge carriers through the weighting potential of 
that electrode. 

The experimental characterisation data sets have been 
used to adjust and optimise the parameters of the code. 
Fig. 38 shows an example comparison between the aver- 
age experimental pulses collected from AGATA crystal 
COOl and the ones calculated with MGS. 

The JASS toolkit [59] provides a Java-based pulse- 
shape simulation using a very similar logic structure as 
MGS shown in Fig. 36. In addition, it allows for a much 
higher granularity in the grid as well as a more pre- 
cise description of the complex geometry of the AGATA 
detectors and the processes at the segment boundaries. 
In connection with an extremely accurate interpolation 
method the description is very precise, especially in the 



critical regions of the crystals. With its automated gen- 
eration of the signal basis for the AGATA detectors an 
excellent agreement between simulated signals and the 
experimentally characterised detector signals has been 
found. 

The ADL package [60-63] allows for the simulation 
of the position-dependent detector response to y-ray in- 
teractions. This library was written in C and it is used 
to create a pre-calculated set of position-sensitive de- 
tector responses. Such sets are used as a lookup table to 
translate online acquired signals into position informa- 
tion. ADL makes use of realistic anisotropic electron- 
and hole-mobility models specially developed for ger- 
manium. The code has a unique capability to simulate 
partially depleted detectors and their capacitances and 
this feature allows for the space-charge reconstruction 
in highly-segmented detectors from capacitance-voltage 
measurements [67, 68]. The code was recently also suc- 
cessfully used to simulate the position-dependent col- 
lection efficiency within the AGATA detectors. Such 
efficiencies can be used to correct for trapping effects 
by exploiting the high position sensitivity of the detec- 
tors [69]. The approach utilises the radial dependence of 
the trapping magnitude which allows a parametric cor- 
rection to be applied. 

A code based on the chain of free open-source soft- 
ware OpenCASCADE [70], gmesh [71], libmesh [72], 
and GSL [73], called AGATAGeFEM [74] has also 
been developed. The open-source software are used to 
model the geometry, to create a mesh adapted as a func- 
tion of the field shape for the solving of the partial differ- 
ential equations, to solve these equations, and finally to 
solve the problem of transporting the charges inside the 
crystal and the formation of the signals. The response 
of the preamplifier and crosstalk are included. A pulse- 
shape database calculated using this code has yet to be 
tested with real data. 

10. Pulse-shape analysis 

The task of pulse-shape analysis (PSA) is to identify 
with high precision the location of the individual inter- 
action points and the corresponding energy deposits of 
a y ray. A y ray will normally have a chain of interac- 
tions in the shell of germanium detectors (e.g. 3 to 4 
at 1.3 Me V). There can be more than one interaction in 
one detector segment and/or the y ray can be scattered 
to another segment of the same crystal or to an adjacent 
detector, or even across the shell. The accuracy of the 
location of these interaction points has to be better than 
5 mm (FWHM). Such a high precision is required in or- 
der to perform the tracking process with high efficiency. 
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Figure 38: (Colour online) A figure showing the 37 charge pulses from the COOl detector for an interaction in a "true coaxial" region of the crystal 
in segment d5. The average experimental pulses are shown in red, the MGS calculated responses are shown in blue and the difference in magnitude 
between the two pulse-shape data sets is shown in black. The difference is shifted down by 0. 1 for clarity. These data were obtained from a 
coincidence scan using a '^^Cs source of 662 keV y rays. The AGATA detector records the 374 keV y rays that have been Compton scattered by 
90°. 



The location and energy determination must be per- 
formed by algorithms that are fast enough for real time 
application, with the computing power available. In- 
deed, the high experimental event rate multiplied by the 
data rate of the digitiser produces an information flux 
far beyond the realistic disk storage capacities in terms 
of volume and recording speed rate. Data from a set of 
37 X 100 ADC samples delivered by the digitiser are re- 
duced through PSA to few parameters of interest. These 
are determined, by comparing the detected pulse shapes 
to a calculated reference basis, in real time. This pro- 
vides the three dimensional interaction position, energy, 
time and a confidence in the quality of the determined 
fit. 

A set of algorithms, such as grid search [57, 75], ge- 
netic algorithms [76], wavelet decomposition, and a ma- 
trix method [77, 78], optimised for diff'erent types of 



events, have been developed. A first operational version 
of the PSA code, a stepwise refined grid-search [57], 
has already been implemented in the final Narval en- 
vironment. Other faster algorithms, which have been 
extensively tested with simulated data and are suited to 
more complex event structures, will be implemented in 
the near future. It is expected that the performance of 
PSA will continuously improve during the project due 
to refined algorithms and increased computing power 
but especially due to improvements in the quality of the 
reference signals. 

The approach based on a genetic algorithm aiming 
for a full decomposition was found to be too slow for 
online processing [76]. Therefore, the collaboration fo- 
cused initially on the grid-search algorithm [57], which 
due to its simplicity is the most robust among the signal 
decomposition codes. 
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The code based on the grid-search algorithm has been 
validated with experimental data obtained in two in- 
beam experiments using techniques based on Doppler 
effects of y rays emitted by nuclei in flight [79, 80]. The 
first experiment was performed at IKP Cologne with 
the prototype ATC detector, which consisted of three 
symmetric HPGe crystals [79]. The second experiment 
was the first in-beam AGATA commissioning experi- 
ment performed at LNL with the first production se- 
ries ATC detector, which consisted of three asymmetric 
HPGe crystals [80]. The measured interaction position 
resolution versus y-ray energy obtained in the latter ex- 
periment is shown in Fig. 39. The grid-search algorithm 
has also been validated in a measurement with the setup 
at LNL using a ^^Na source [81]. 
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Figure 39: Interaction position resolution (FWHM) versus y- 
ray energy Ey extracted from the data obtained in tlie first AGATA 
in-beam commissioning experiment based on an inverse kinematics 
fusion-evaporation reaction with a 64MeV ^^"Si beam and a thin '^C 
target. The error bars include both statistical and systematic errors. 
See [80] for further details. 

It was shown that in the case of single-hit events the 
grid-search algorithm already gives a sufficiently good 
Doppler correction. However, it is not well suited to 
more complex events, the reason being that in the case 
of multiple interactions in a crystal, or especially in 
a segment, the analysis effort strongly increases with 
the number of degrees of freedom. However, the grid- 
search algoritm has been used successfully with in- 
beam data to handle multiple interactions in the same 
crystal as along as the segments with transient signal do 
not overlap. 



Several fitting procedures as mentioned above have 
been developed in order to deal with these types of 
events. The currently favoured two methods are, 
the fully informed particle swarm [75] or the SVD 
(singular-value decomposition) matrix method [77, 78, 
82, 83]. The former method is faster but the latter gives 
a better precision on the location of the hits. All these 
new algorithms still have to be validated with AGATA 
data. Given the absence of a one-for-all algorithm a 
dispatcher code, choosing the most adequate PSA algo- 
rithm for a given event, will also be implemented soon. 

Table 10 gives an overview on the performance of 
different algorithms derived from tests using simulated 
data. The CPU performance is quoted relative to the 
grid-search algorithm (3 ms per event) for single inter- 
actions and the particle-swarm algorithm (2 ms to 5 ms 
per event) for multiple interactions. 

After the energy calibration, which is done by the 
front-end electronics, the influence of noise, pedestal, 
time jitter and cross talk needs to be understood in or- 
der for the experimental signals to be successfully com- 
pared with the basis signals 

Concerning the time jitter, different time-shift deter- 
mination methods have been considered: supplemen- 
tation of the basis with time-shifted signals. Neural 
Network determination [75], Kolmogorov-Smirnov de- 
termination [86], Taylor expansion and the substitu- 
tion of residue minimisation by chi-square minimisa- 
tion [87, 88]. The first method has the drawback of 
requiring larger signal bases, the latter method reduces 
the influence of the time jitter but does not measure 
it. Using simulated data, the Neural Network algo- 
rithm has shown to be very fast and provides a reso- 
lution better than 2 ns at computing times in the order of 
10 [.IS. Taking into account only the core signal, this al- 
gorithm seems to be very robust and hardly affected by 
noise and cross talk discussed below. The Kolmogorov- 
Smirnov and Taylor expansion methods are purely alge- 
braic. They also give excellent results in terms of com- 
puting speed, robustness and precision. 

The cross-talk effect is present in any segmented de- 
tector. It induces energy shifts and decreases the hit- 
location precision as it mixes the transient signals. Tran- 
sient signals are those induced on adjacent electrodes 
to the primary interaction due to the drift of the charge 
carriers inside the germanium crystal. Cross-talk con- 
tributions can also appear between segments of different 
detectors. Only negligible cross-talk components were 
observed between a firing detector and the core signal 
of a neighbouring detector [22, 62]. The effects of cross 
talk are described in details in section 3. Two quite sim- 
ilar correction methods have been used. The first one 
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Table 4: Overview of the performance of different PSA algorithms derived from tests using simulated data. The performance is quoted relative to 
the grid-search algorithm. The position resolution is quoted in as FWHM values. 





Single interaction 


Multiple interactions 


Algorithm 


CPU time 


Resolution 


CPU time 


Resolution 




[norm/GS] 


[mm] 


[norm/PS] 


[mm] 


Grid search (GS) [57] 


1 


2 






Extensive GS [75] 


2.7 


1 


6-104 


4 


Particle swarm (PS) [75] 


0.1 


2 


1 


4 


Matrix method [77, 78] 


6.7 


2.4 


10 


10 


Genetic algorithm [76] 


330 


1.9 


2x 10^ 


8.1 


Binary search 


0.06 


1 


Not adapted 


Not adapted 


Recursive subtraction [84, 85] 


Not evaluated 


3" 


Not evaluated 


5^ 


Neural network [75] 


2 


1.5 


Not adapted 


Not adapted 


Wavelets 


Not evaluated 


2.3 


Not adapted 


Not adapted 



^ Only neighbouring segment hits. 
Only radial coordinate. 



includes the cross-talk effect into the basis signal using 
a cross-talk matrix. The second corrects the measured 
signal by the cross-talk effect. These corrections per- 
mit the lowering of the FWHM error on the position by 
about 1 mm. 

Another difficulty in the determination of the precise 
interaction position arises from the fact that the germa- 
nium crystal may be slightly displaced or tilted inside 
its capsule. It has been shown that this problem may be 
addressed prior to the experiment using the scanning ta- 
bles or within the final setup using a radioactive source 
at fixed points and the y-imaging method [89]. 

In order to facilitate the maintenance and develop- 
ment of algorithms used for AGATA two levels of ab- 
straction are used for the PSA. The first is the ADF li- 
brary [44] that is responsible for coding and decoding 
data in the data flow. The second layer of abstraction is 
a set of C++ classes, which provide, via inheritance, sim- 
ple means of "attaching" a PSA algorithm to the Narval 
DAQ. This system has been successfully tested with the 
grid-search algorithm. Due to the fact that no cross cor- 
relations between the signals of different AGATA crys- 
tals has been found so far, all the PSA relevant algo- 
rithms operate on the signals of each individual crystal 
in parallel on different machines. This allows the use of 
several computers in parallel to analyse different events. 
Currently, it is foreseen to use 2 processors with 4 cores 
each to process the data of one crystal which will im- 
prove the PSA performance by a factor of 2 to 4. 

Finally the events have to be reassembled according 
to their timestamps and a tracking algorithm is apphed 



in order to disentangle the coincident interaction points 
and to determine the total energy and the emission di- 
rection of those y rays that have been fully absorbed in 
the germanium shell. Absolute positions of the individ- 
ual crystals, tilting angles and target position corrections 
enter at this stage. 

11. Gamma-ray tracking 

The aim of tracking algorithms is to reconstruct the 
trajectories of the incident photons in order to determine 
their energy and direction. To do this, the algorithms 
must disentangle the interaction points identified in the 
detectors and establish the proper sequences of interac- 
tion points. Tracking algorithms can be divided into two 
classes: algorithms based on back tracking [90] and al- 
gorithms based on clustering and forward tracking [91]. 
Both are related to the particular properties of the in- 
teraction of photons with matter. In forward-tracking 
algorithms, the first step of the procedure is to group in- 
teraction points into clusters in (6, </>) space. Back track- 
ing, on the other hand, is based on the observation that 
the final photoelectric interaction after a y-ray scattering 
sequence usually falls into a narrow energy band. Start- 
ing from interaction points in that energy range the al- 
gorithm tracks back towards the original emission point 
using the physical characteristics of the interaction pro- 
cess and selects the most probable interaction scheme. 

For photon energies of interest (from about lOkeV 
to about 20 Me V), the main physical processes that oc- 
cur when a photon interacts in germanium are Comp- 
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ton scattering, Rayleigh scattering, pair creation and the 
photoelectric interaction. Since Compton scattering is 
the dominant process between 150keV and lOMeV, all 
tracking algorithms are based on the properties of this 
interaction. 

Within the AGATA project, the development and op- 
timisation of tracking algorithms have been performed. 
Tracking has also been used to evaluate the performance 
of PSA. Another part of the work related to tracking has 
been the use of tracking techniques to investigate the 
background-suppression capabilities from sources of ra- 
diation that do not originate from the target. 

11.1. Development of tracking algorithms 

The development of tracking algorithms and the im- 
provement of existing ones is an ongoing process. Re- 
alistic simulated data sets, produced with the Geant4 
AGATA code [14], have been and are systematically 
used to test, compare and improve the performance of 
tracking algorithms. This code treats the physics of the 
energy deposition correctly by generating and propagat- 
ing all the possible secondary particles and by taking 
into account the momentum profile of electrons in ger- 
manium [92]. In order to produce more realistic data, 
the interaction points separated by less than c/jes = 5 mm 
are packed together, the positions of the interaction 
points are randomly shifted in all direction {x, y and z) 
by sampling a Gaussian energy-dependent uncertainty 
distribution and an energy threshold of 5 keV is applied 
to all the interaction points within each segment. 

It has been shown that forward-tracking algorithms 
are more efficient than back-tracking ones [92]. This 
is because back-tracking algorithms rely on the identi- 
fication of the last (photoelectric) interaction point in a 
scattering sequence, which generally loses its original 
characteristics after PSA: it is poorly localised and/or 
packed with other interaction points. However, ulti- 
mately the optimal y-ray tracking algorithms may re- 
quire the use of a combined approach where for exam- 
ple back tracking is used to recover events missed by the 
clustering techniques. 

77.2. Clustering techniques 

In the Mars Gamma-Ray Tracking (MGT) [93] and 
the Orsay Forward Tracking (OFT) [92] codes, both 
based on the forward-tracking technique, points are 
grouped into clusters according to their relative angu- 
lar separation. In order to increase the tracking perfor- 
mance for high photon multiplicities without decreasing 
the performance at lower multiplicities, the maximum 
allowed separation angle is set to depend on the number 
of interaction points in the event. 



A clustering technique based on the Fuzzy C-Means 
algorithm has been developed [94]. After running this 
algorithm, the optimal positions of the centres of the 
clusters and the degree of membership of each point to 
each of the clusters are obtained. The next step in the 
code is the "defuzzification", making possible the use 
of standard validation procedures. Another clustering 
method, called Deterministic Annealing Filter (DAF), 
has also been developed [95]. In this method, the clus- 
tering problem is to perform the minimisation of the dis- 
tortion criterion, which corresponds to the free energy 
related to the annealing process in statistical physics. 

The detection efficiency and peak-to-total obtained 
by the above-mentioned algorithms for 1 MeV photons 
emitted at the centre of AGATA are summarised in Ta- 
ble 5. 

Table 5: Simulated full-energy efficiency and peak-to-total (P/T) of 
AGATA for cascades of 30 (and 1)1 MeV photons. In all the cases, 
the data is packed and smeared in the standard way and a 5 keV energy 
threshold has been applied. 



Algorithm 


Efficiency 


P/T 




[%] 


[%] 


MGT 


28(43) 


49(58) 


OFT 


24(37) 


54(68) 


Fuzzy C-Means 


27 


46 


DAF 


26(36) 


47(69) 



77.5. Effective distances 

In all tracking philosophies, the photon trajectories 
are tracked with the help of the Compton scattering for- 
mula and with the ranges and cross sections of the pho- 
toelectric, Compton and pair-production interactions. 
To compute the appropriate ranges of photons in ger- 
manium, the effective distances in germanium between 
interaction points as well as between the interaction 
points and the source need to be computed. The prob- 
lem can be solved analytically if the detector geometry 
is approximated to a An germanium shell. This approx- 
imation was validated with the Geant4 simulation code 
[14], in which the exact geometry of AGATA is defined, 
i.e. the shape of the crystals, encapsulation, cryostats, 
the empty spaces and the distances between all these el- 
ements. 

77.4. PSA and front-end electronics 

It has recently become clear that PSA algorithms 
have difficulties in dealing with events in which there 
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is more than one interaction point per segment. This 
is because the parameter space to search becomes very 
large and because the properties of the response func- 
tion of the detector can lead to more than one solution. 
Although algorithms to determine the number of inter- 
action points in each segment have been developed and 
have shown encouraging results in the coaxial part of 
the crystal [85], the effect of systematically assuming 
only one interaction point per segment has been imple- 
mented. No drastic drop in performance is observed, 
not even when the recoil velocity is 50 % of the speed 
of light, which yields a higher concentration of interac- 
tion points because of the Lorentz boost. This result is 
understood since the average of the distribution of the 
number of interaction points per segment is close to 1 
and varies only slightly with incident photon energy. 

Due to electronic noise, signals corresponding to a 
deposited energy below a certain value will not be de- 
tected in the segments of AGATA. The effect of this en- 
ergy threshold on the tracking efficiency has been in- 
vestigated. The loss in full-energy efficiency with re- 
spect to the case of a 1 keV threshold is found to be 
1% (5keV threshold), 8% (lOkeV threshold), 13% 
(20keV threshold) and 29% (50keV threshold). The 
data from AGATA has shown that most of the full- 
energy efficiency can be recovered by using the core- 
energy information. 

11.5. Neutrons 

The influence of the elastic and inelastic interactions 
of neutrons on tracking performance has been found to 
be significant for low to medium y-ray multiplicities es- 
pecially on the P/T ratio [96]. As far as pulse shapes are 
concerned, it has been shown that neutrons and y rays 
yield similar signals in nonsegmented closed-end coax- 
ial and planar «-type HPGe crystals, even though the in- 
teraction mechanisms are different [97]. Distinguishing 
neutrons and y rays by time of flight methods requires 
a timing resolution better than about 5 ns, which will 
be hard to achieve in the case of low-amplitude signals. 
Fingerprints of the neutron interaction points have been 
found, which may be used to reduce the background due 
to neutrons without too much loss in full-energy effi- 
ciency [98]. 

11.6. Sea tte ring ma te ria Is 

The effect of passive materials in the inner space of 
the array, especially in connection with the used of com- 
plementary instrumentation, has been thoroughly dis- 
cussed in Ref [14]. The conclusion presented in this 



work is that at medium or high y-ray energies, the pres- 
ence of scattering materials does not affect the perfor- 
mance of a tracking array more than in the case of a 
conventional array. However, at low energy Rayleigh 
scattering changes the direction of the incident photon, 
which directly affects the tracking process. 

77.7. Imaging 

The performance of PSA is generally tested by check- 
ing how well experimental spectra can be Doppler- 
corrected since this procedure directly depends on the 
attainable precision in locating each interaction point. 
Another method has recently been devised. It relies 
on Compton imaging, which does not need any beam 
nor complex experimental setup. Instead it is assumed 
that the quality of the image reconstruction obtained 
through the knowledge of the Compton scattering se- 
quence in the detector can provide information on the 
interaction position resolution. The position and en- 
ergy of the interaction points were extracted using the 
grid-search algorithm [57] assuming only 1 interaction 
point per segment. Since the incident energy is known, 
the first scattering angle can be extracted by selecting 
events according to their total energy and number of in- 
teraction points. This angle defines a scattering cone in 
3D space. If more than one event is analysed, all the 
cones should overlap in a single point corresponding to 
the source of the y rays. Placing the source far away 
from the detector reduces the problem to the forming 
of a 2D image on the surface of a sphere, which can 
be represented in a {9, <p) plane. Comparison of the ex- 
perimental and simulated and (p image profiles of a 
^"Co source has yielded an interaction position resolu- 
tion of about FWHM - 5 mm [99], which is compara- 
ble with the result obtained from in-beam experiments 
[79, 80, 100]. Attempts, using simple algorithms, to 
evaluate the capability of the AGATA detectors to reject 
y rays originating from locations that are different from 
the target position have shown partial success [101]. 

11.8. Integration into the DAQ 

Tracking is performed on 1 or 2 nodes of the PC farm 
of the AGATA DAQ, which runs the Narval system. 
The integration of the tracking and other algorithms into 
Narval were done by using the ADF library [44]. The 
goal is to make the data format transparent to the algo- 
rithms, i.e. ADF allows each algorithm to access the 
data it needs without knowing the structure of the data 
flow. It also provides a virtual Narval environment in 
which to test and debug algorithm codes in stand-alone 
mode. OFT and MGT are currently the only tracking 
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algorithms, which have been integrated into the Narval 
architecture. The online performance of OFT has been 
optimised and meets the specifications. 

77.9. Tracked Co spectrum 

The results of a measurement with a ''"Co source with 
three ATC detectors are shown in Fig. 40. The source 
was placed in the focal point of AGATA at a distance of 
235 mm from the front face of the crystals. The black 
histogram shows a "raw" crystal spectrum obtained by 
summing all the individual gain-matched energy depo- 
sitions in the crystal, corrected for segment threshold ef- 
fects with the core-energy signal. The blue histogram is 
a tracked spectrum obtained by OFT and by excluding 
all single interaction points. The PSA was performed 
using the grid-search algorithm (see section 10) with the 
ADL pulse-shape database (see section 9). The P/T of 
the raw spectrum is 16.8 % using a low-energy cut-ofF 
of OkeV. The P/T of the tracked spectrum is 53.9% 
and 54.6 % using OkeV and 200 keV cut-off values, re- 
spectively. The tracking efficiency, defined as the ra- 
tio of the number of counts in the full-energy peak in 
the tracked spectrum to the number of counts in the raw 
crystal spectrum, is 84.5 %. 
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Figure 40: (Colour online) ^^Co spectrum measured with three ATC 
detectors. Black histogram (greater number of counts): obtained by 
using the gain matched core signal of the crystals. Blue histogram 
(fewer number of counts): tracked spectrum excluding single interac- 
tion points. See text for further details. 



12. Data analysis 

As described in section 8, at the end of the AGATA 
DAQ chain, coincidences between tracked y rays and 
ancillaries are provided. Data analysis needs to be per- 
formed at several stages, including replaying the full ex- 
periment from the raw traces. To help in these tasks a 



package called GammaWare (GW) [47] has been devel- 
oped. 

The GW package is divided into several sub-projects; 
Core, Physics, GEM (Gamma-ray Event Monte Carlo), 
Tools, ADF, ADFE (Extended ADF library). A dy- 
namic library, written in C++, is associated to each spe- 
cific part. All libraries, except ADF, are based on the 
ROOT framework [102]. As a collaborative software, 
GW is available through the open-source version con- 
trol system subversion and uses a bugtracker system. 
User's documentation and web documentation, made 
by doxygen, can be found at the AGATA data-analysis 
web site [47]. 

The GW package adds facilities to ROOT that are 
specific to analysis of y-ray spectroscopy data. Only the 
most significant add-ons are described here. The Core 
library contains facilities common to all sub-projects. 
The Physics library defines C++ objects specific to y- 
ray spectroscopy analysis such as level schemes. As 
any ROOT object, a level scheme can be displayed and 
saved in ROOT files. They can be built graphically or 
imported from ENSDF [103] and RadWare [104] files. 
As an example, using the graphical level scheme dis- 
play, gates are selected and then applied on correlated 
spaces, which are (virtual C++) interfaces to any kind 
of system storing y-ray coincidences (y"). So far y-y 
(y^) like matrices are implemented but the system is de- 
signed to be extendable to higher coincidence orders. 
From level schemes, cascades of discrete y rays can be 
randomly generated and these cascades can be fed into 
Geant4 simulations: the GEM sub-project is devoted to 
physics generators. In the Tools library many useful fa- 
cilities can be found. It is the place where links to other 
data analysis frameworks are implemented, allowing for 
instance to ease exchanges of ID and 2D spectra. 

The GW package includes the ADF library (see sec- 
tion 8), which can be used also in stand-alone mode. 
Such a library contains a complete virtual interface to 
any actor processing the AGATA data flow. It defines 
the format of the data that are exchanged between the 
AGATA algorithms: they are encapsulated in frames 
with a trigger mechanism to deliver to a particular al- 
gorithm only the specified frames. 

All the C++ actors processing the data flow are 
grouped in a package consisting of several dynamic li- 
braries that are loaded into Narval. The actors can 
also be used in any other framework, in particular 
ROOT/GW, as well as linked to build a standard ex- 
ecutable program. Some actors can process the data 
flow and moreover perform first stage analyses (count- 
ing rates, histograms, even traces saved in ROOT or 
ASCII files) to check online that the system is running 
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Figure 4 1 : (Colour online) Online tracking Watcher display showing 
typical Compton scattering events in the 15 Ge crystals of AGATA, as 
provided by the tracking algorithm. Only the first and second interac- 
tion points are displayed. 



properly. 

In the GW package, the ADF library is extended 
through the ADFE sub-project bringing the AGATA 
like data to the ROOT/GW environment. In particular 
Watchers are defined there. They are light tasks dedi- 
cated to specific frames (output of an algorithm) running 
independently and easy to plug whatever the data flow 
structure is. A typical goal of the Watchers is to build 
histograms/graphs, to calibrate events from ancillary de- 
tectors. It is also used to save data in ROOT trees and 
to display events in 3D geometries via the ROOT/EVE 
facility. This is illustrated in Fig. 41 in which Comp- 
ton events as provided on-line by the tracking algorithm 
can be visualised. The spy mechanism proposed by Nar- 
val (section 8) allows the same Watchers to be used on- 
line/offline and allows analysis to be performed at all 
stages of the data flow: signals (core and segments), hits 
distributions, raw data from the AGAVA card, counting 
rates, etc. 

The virtualisation provided by the ADF library al- 
lows for the development of emulators, i.e. actors con- 
nected together within a common environment in vari- 
ous topologies as it is done in Narval [45]. Contrary to 
Narval, emulators only run so far on single computers 
and are not efficient in case of complex situations de- 
spite the fact that some of them run in multi-threading 
mode. However, as they are written in C/C++, it makes 
the development, testing, and debugging phases of the 
AGATA data-processing algorithms easier Most of the 
actors are shared libraries that can be plugged in Narval 



or in emulators. For efficiency reasons the event builder 
and the merger running online are written in Ada so they 
cannot be used as such for offline purposes. To over- 
come this problem, specific C++ actors have been devel- 
oped to reproduce event builders and mergers as they 
act online. Emulators provide also elegant solutions to 
build almost any kind of complex analysis chain. Sim- 
ple ones are used in the GW package to launch Watch- 
ers. More complete ones have been developed as for 
example the emulator of the whole topology as defined 
for the actual experiments. 

72.7. Data processing on the Grid 

The AGATA collaboration has adopted the Grid tech- 
nologies [105] and has created its own virtual organi- 
sation (VO), namely vo.agata.org. The Grid computing 
resources (tape, disk and CPU) are presently provided 
by Lyon, Strasbourg, Paris and Valencia and are shared 
and accessed by the VO members under agreed com- 
mon policies. 

As mentioned in section 8, the raw data containing 
the preamplifier traces are recorded on a Grid tape- 
storage system at Tierl sites for future off-line re- 
play using, for instance, more sophisticated PSA and 
tracking algorithms than those implemented to run on- 
line. About 150 TB of data have been produced during 
the commissioning and the first experiments performed 
with up to 15 crystals during the period from March 
2009 to September 20 11, with an average value of 8 TB 
of data collected per experiment. 

Besides the storage space question, partly discussed 
in [106], the transfer of the raw data to local institutes 
for analysis, as well as their processing, are also critical 
issues as they are very time and resource consuming. 
As the raw data have been stored on the Grid, then it is 
advisable to use also the Grid computing power (CPU) 
to replay the data, by emulating all or part of the DAQ 
system processing stages. 

It has been demonstrated that replaying raw data on 
the Grid using an emulator is feasible without the need 
of changing the software to run on the Grid [107]. This 
is the case as long as copies of the input data files are 
previously downloaded into the node where the emula- 
tor is running, so the data are accessed locally. However, 
in the framework of designing a Grid computing model 
for data management and data processing for AGATA 
[108], more developments have been recently tested in 
order to access directly the data on the Grid storage 
[109]. For this, the emulator software has been modi- 
fied by including the necessary functions from the Grid 
File Access Library (GFAL) in order to interact with the 
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Grid storage resource managers, through the appropri- 
ate protocols. 

In the following, an example is presented in order 
to illustrate the performances in execution times for re- 
playing raw data on the Grid. Three tasks of 49, 70 and 
140 jobs, respectively, have been run on the IFIC/GRID- 
CSIC e-science infrastructure, each job executing an in- 
stance of the emulator that processes PSA and tracking 
on a fraction of the considered raw data set. It is worth 
noting that this infrastructure provides a distributed file 
system, namely Lustre, which allows direct access to the 
data located at the Grid storage of the site, without us- 
ing GFAL. A total of respectively 2, 3 and 6 TB of com- 
missioning data (147, 210 and 420 data files of 14 GB 
each, about 10** events) have been processed. Fig. 42 
shows the evolution with time of the execution of the 
jobs for the above defined tasks. The average values 
of time measured for the completeness of the tasks, ob- 
tained from re-submitting them to Grid several times, 
are 3, 7 and 15 hours to process 2, 3 and 6 TB data, 
respectively. These results indicate that 8 TB of data, 
produced in a typical experiment, can be processed in 
less than one day. 
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Figure 42: Completed jobs of data processing for 3 tasks run on the 
Grid as a function of time. Each job runs an instance of the emulator 
to replay part of the AGATA raw data set of interest. 



13. Summary and outlook 

The realisation of the AGATA spectrometer is a result 
of many technological advances. These range from ad- 
vances in Ge detector technology, digital DAQ systems, 
signal decomposition and y-ray interaction reconstruc- 
tion, and in many areas of the infrastructure needed to 
support and operate such a complex device. 



The AGATA spectrometer is now fully operational 
in its first physics campaign at INFN LNL in Legnaro, 
Italy, utilizing the wide range of stable beams available. 
A view of AGATA at the target position of the PRISMA 
spectrometer is shown in Fig. 43. AGATA is designed to 
be a peripatetic instrument and will move between ma- 
jor laboratories in Europe to take advantage of the range 
of different beams and equipment at each laboratory and 
of the resulting scientific opportunities. AGATA will 
be operated in a series of campaigns, the one first at 
LNL, and subsequently at the GSI facility in Germany 
and the GANIL laboratory in France. At GSI, AGATA 
will be used at the exit of the Fragment Separator (FRS) 
to study very exotic nuclei produced following high- 
energy fragmentation and secondary Coulomb excita- 
tion. At GANIL it will use the wide range of radioactive 
ions from the coupled cyclotrons and SPIRAL. During 
these first three physics campaigns the array will con- 
tinually increase its efficiency as more detector systems 
are added. The first stage is to build up the system to 
60 detector crystals, and then proceed towards the full 
implementation of the An AGATA. Subsequent physics 
campaigns will take advantage of the new radioactive 
beams available as facilities such as FAIR, SPIRAL2, 
SPES and HIE-ISOLDE come onUne. 

AGATA will have an enormous impact on nuclear 
physics research in particular the exploration of nuclear 
structure at the extremes of isospin, mass, angular mo- 
mentum, excitation energy, and temperature. This rad- 
ically new device will constitute a dramatic advance in 
y-ray detection sensitivity that will enable the discov- 
ery of new phenomena in nuclei, which are only pop- 
ulated in a tiny fraction of the total reaction cross sec- 
tion or that are only produced with rates of the order of 
a few per second or less. The unprecedented angular 
resolution afforded by its position sensitivity will facili- 
tate high-resolution spectroscopy with fast and ultra-fast 
fragmented beams giving access to the detailed structure 
of the most exotic nuclei that can be reached. In addi- 
tion, the capability to operate at much higher event rates 
will allow the array to be operated for reactions with 
intense y-ray backgrounds. 

The instrumentation and technical advances driven by 
this work, and the knowledge gained by those involved, 
is also important in a wide range of applications. These 
advances have potential impact in areas such as medi- 
cal imaging systems, homeland security, environmental 
monitoring and the nuclear industry. 

In addition to the technical advances, AGATA rep- 
resents a tremendous human achievement in the suc- 
cessful collaboration of several hundred personnel in 12 
countries and over 40 laboratories and institutes across 
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Figure 43: (Colour online) The AGATA and PRISMA spectrometers at LNL. The beam from the accelerators enters from the left. A part of the 
beam line, close to AGATA and half the scattering chamber has been removed. The quadrupole and dipole magnets of PRISMA are seen on the 
right-hand side of the photograph. 



Europe. The collaboration is now excited with the 
prospect of using this spectrometer and to capitalise 
on its discovery potential for the understanding of the 
atomic nucleus. 
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